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ABSTRACT 

Wc present a new distance estimation method for dust-continuum-identificd molecular 
cloud clumps. Recent (sub-)millimeter Galactic plane surveys have cataloged tens of thou- 
sands of these objects, plausible precursors to stellar clusters, but detailed study of their 
physical properties requires robust distance determinations. We derive Bayesian distance 
probability density functions (DPDFs) for 770 objects from the Bolocam Galactic Plane 
Survey in the Galactic longitude range 7?5 < ^ < 65°. The DPDF formalism is based on 
kinematic distances, and uses any number of external data sets to place prior distance prob- 
abilities to resolve the kinematic distance ambiguity (KDA) for objects in the inner Galaxy. 
We present here priors related to the mid-infrared absorption of dust in dense molecular 
regions and the distribution of molecular gas in the Galactic disk. By assuming a numer- 
ical model of Galactic mid-infrared emission and simple radiative transfer, we match the 
morphology of (sub-) millimeter thermal dust emission with mid-infrared absorption to com- 
pute a prior DPDF for distance discrimination. Selecting objects first from (sub-)millimeter 
source catalogs avoids a bias towards the darkest infrared dark clouds (IRDCs) and ex- 
tends the range of heliocentric distance probed by mid-infrared extinction and includes 
lower-contrast sources. We derive well-constrained KDA resolutions for 618 molecular cloud 
clumps, with approximately 15% placed at or beyond the tangent distance. Objects with 
mid-infrared contrast sufhcient to be cataloged as IRDCs are generally placed at the near 
kinematic distance. Distance comparisons with Galactic Ring Survey KDA resolutions yield 
a 92% agreement. A faee-on view of the Milky Way using resolved distances reveals sections 
of the Sagittarius and Scutum-Centaurus Arms. This KDA-resolution method for large cat- 
alogs of sources through the combination of (sub-) millimeter and mid-infrared observations 
of molecular cloud clumps is generally applicable to other dust-continuum Galactic plane 
surveys. 

Subject headings: Galaxy: kinematics and dynamics, structure - infrared: ISM ~ ISM: 
clouds, dust - stars: formation 



I. INTRODUCTION 

Recent (sub-)milli meter surveys of the Galactic 
plane (ATLASGAL ISchuller et aH [200l Hi -GAL 
IMolinari et all [2010: BGPS. lAguirre et al.l l20Tll) 
have detected tens of thousands of molecular cloud 
cores and clumps in thermal dust emission. As 
plausible precursors to stellar clusters, OB associ- 
ations, or smaller stellar groups, molecular cloud 
clumps can yield clues about the for mation of mas- 
sive stars (jMcKee fc Ostrikerll2007D . The masses 
and temperature profiles of these objects are key to 
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unraveling this process. Rece nt work has sought 
to measure thes e quantities (jRusseil et al.l 120111 : 
lEden et al]l2012[ ). but a robust and comprehensive 
tally does not yet exist. 

Derivation of masses for molecular cloud clumps 
from dust continuum data requires an estimate of 
the heliocentric distance to each object and the tem- 
perature of the emitting dust. Analysis of Her- 
schel Hi-GAL data is beginning to yield te mper- 
ature maps of the Galact ic plane l^eretto et al.l 
I2010t lBattersbveraIll2011[ ). While a detailed un- 
derstanding of the interplay between dust temper- 
ature and the environment and evolution of molec- 
ular cloud clumps is important, variations in the 
assumed dust temperature by a factor of two only 
produce a factor of a few difference in the mass de- 
rived from (sub-)millimeter observations. In con- 
trast, the derived mass of a molecular cloud clump 
is proportional to the square of its heliocentric 
distance; accurate distance estimates play a far 
larger role in the mass calculation. Recent studies 
of isolated regions, with well-det ermined distances , 
such as Perseus and Ophiuchus (iRidge et edl 120061 : 
lEnoch et all 120061: iRosolowskv et al.l 12008( 1 have 
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unveiled ma ny properties of mo l ecular cloud cores i n 
recent years (|Enoch et al . 11200 TtlSchnee et al.ll2010f ). 
To gain similar insight into the larger molecular 
cloud clumps seen spread throughout the Galactic 
plane, a robust method for distance determinations 
for large data sets is required, because the distances 
to most clumps are subject to the kinematic dis- 
tance ambiguity (KDA). 

The most straightforward method for estimat- 
ing the heliocentric distance (dg) to a molecu- 
lar cloud clump is to project its observed linc-of- 
sight velocity (wlsr ) , derived from molecular line 
Doppler shifts, onto a Galactic rotation curve. 
These kinematic distances are generally unique for 
the outer Galaxy, but inner Galaxy sources are 
subject to the KDA, a projection effect of the or- 
bital motion for objects within the Solar Circle 
(i?o)- A line of sight intersecting a circular or- 
bit at Galactocentric radius i?gai < Rq crosses 
that orbit twice, each with different spatial ve- 
locities but both with the same Wlsr- Various 
techniques have been suggeste d for resolving the 
KDA ( 21-cm H I absorpti o n: [ Anderson fc Baiiial 
120091 IRoman-Duval et all 120091: the presence 
of mid-infrared dark clouds: iRathborne et alj 



200l iPeretto fc Fulled [200l H2CO absorption 



Sewilo et al. _ 20041: and near-infrared extinction: 



Marshah et al.ll2009llFoster et al.ll2012[) : this paper 



presents a method based on comparing mid-infrared 
extinction with (sub-)millimeter emission. 

Appearing as dark absorption features against 
a bright mid-infrared background, infrared dark 
clouds (IRDCs) offer a practicable means for re- 
solving the KDA. IRDCs are most striking against 
the broad, diff use Gala ctic emission nea r A = 8 /im 
(jPerault et all [1996: S imon et al.ll2006[ ). although 
they may be detected in absorption against back- 
ground stars at o ther infrared wavelengths (c./. 
[Foster et al.l IMl) . Studies of IRDCs at (sub- 
)millimeter wavelengths reveal that they are dense 
molecular cloud clumps (iJohnst one et al.l 12003 



IRathborne 'eral]\2UU^ IBattCTsbv et al.ll20ia Mm . 
As extinction features, IRDCs must lie in front of 
enough mid-infrared emission to be visible. It is 
possible to a priori assign the near kinematic dis- 
tance for the darkest cloud s {e.g. iButler fc TanI 
2009" 'Pcretto fc Fullei1l2009( ). but recent work by 



Battcrsby ct al] (|2011[ ) has shown that molecular 
cloud clumps may be visible as slight intensity 
decrements in the mid-infrared at the far kinematic 
distance despite not being dark enough to be cat- 
aloged as IRDCs. To encompass this second set 
of objects, we classify all dust-continuum-identified 
molecular cloud clumps with mid-infrared intensity 
decrements of any amount as Eight-Micron Absorp- 
tion Features (EMAFs), whether catalogued as ei- 
ther an IRDC or not. These constitute a general- 
ized collection of cold molecular cloud clumps iden- 
tified first by dust-continuum emission and then 
checked for infrared absorption. The EMAF defi- 
nition excludes objects extensively undergoing the 
later stages of star formation or that are exposed 



to strong ultraviolet radiation, as both processes 
excite PAH emission near A = 8 yLtm, rendering 
invisible any absorption. Investigating the mid- 
infrared properties of molecular cloud clumps based 
on this classification avoids a bias toward the dark- 
est, nearby IRDCs. 

This paper presents a quantitative distance esti- 
mation technique for molecular cloud clumps based 
on Bayes' Theorem. A distance probability den- 
sity function (DPDF) is computed using a distance 
likelihood derived from kinematic information (ob- 
served Wlsr) and prior probabilities, based on an- 
cillary data sets, that are applied in an effort to 
resolve the KDA. We present here two such priors. 
The first involves the comparison between observed 
mid-infrared absorption and millimeter emission of 
individual molecular cloud clumps, and the second 
is based on the Galactic-scale distribution of molec- 
ular gas. In addition to those described here, any 
number of additional priors may be applied to con- 
strain the distance estimate. 

The apparent optical depth of an EMAF calcu- 
lated naively from mid-infrared images is likely less 
than the true value due to diffuse 8-/im emission 
lying between the cloud and the observer. By pa- 
rameterizing the amount of total mid-infrared emis- 
sion along a line of sight lying in front of a molecu- 
lar cloud clump as the "foreground fraction" (/fore), 
simple radiative transfer arguments may be used to 
derive the true optical depth. The recent numerical 
Galactic infrared emission model of lRobitaille et ahl 
()2012|) offers an estimate of /fore as a function of 
dg in the Galactic plane. The maximum likeli- 
hood distance to a molecular cloud clump may be 
derived by comparing the optical depth calculated 
from (sub-)millimeter thermal dust continuum data 
with the absorption optical depth derived from the 
mid-infrared images and ffoi-c{d„). This comparison 
generates a DPDF that takes into account Galactic- 
scale conditions along a given line of sight, including 
spiral structure. A DPDF derived in this manner 
contrasts the widely-used "step-function" method 
whereby a molecular cloud clump is automatically 
assigned the near kinematic distance upon associa- 
tion with a catalogued IRDC. 

The methodology presented here is valid only for 
molecular cloud clumps that exhibit mid-infrared 
absorption, and therefore is but one means for 
distance discrimination for large catalogs of dust- 
continuum-identificd objects. We present an au- 
tomated means for deriving Bayesian DPDFs for 
mid-infrared dark molecular cloud clumps detected 
by the Bolocam Galactic Plane Survey, but this 
method is applicable to all (sub-) millimeter Galac- 
tic plane surveys. 

This paper is organized as follows. Section [2] de- 
scribes the data sets used. The DPDF formalism is 
described in Section [3l Section |4] outlines the gener- 
ation of prior DPDFs for EMAFs. Results from the 
Bayesian DPDFs are presented in Section [S] Impli- 
cations of this work are discussed in Section [6l and 
conclusions are presented in Section [T] 
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2. DATA SETS 
2.1. The Bolocam Galactic Plane Survey 

The Bolocam Galactic Plane Survey (BGPS; 
lAguirreet aL|[20ll1: iGinsburg et"al|[20l3l ) is a A = 
1.1 mm continuum survey covering 170 deg^ at 33" 
resolution. The BGPS was observed with the Bolo- 
cam instrument at the Caltech Submillimeter Ob- 
servatory (CSO) on Mauna Kea. It is one of the 
first large-scale blind surveys of the Galactic Plane 
in this region of the spectrum, covering —10° < £ < 
90° with at least |6| < 0?5, plus selected regions in 
the outer Galaxy. For a map of BGPS VI. cov- 
erage and details about obser vation methods and 
the data reduction pipeline, see lAguirre et al.l (|2011L 
hereafter All). 

From the BGPS Vl.O images, 8,358 millimeter 
dust-continuum sources were identified using a cus- 
tom extraction pipeline. The BGPS catalog (Bolo- 
cat) contains source positions, sizes, and flux den- 
sities extracte d in various apertures, among other 
quantities (see iRosolowskv et al.ll20ldl for complete 
details). BGPS Vl.O pipeline products, including 
image mosaics and the catalog, are publicly avail- 
ablcO. For this work, we utilized the flux densities 
measured in a 40" top-hat aperture, which has the 
same sohd angle as the BGPS 33" FWHM Gaus- 
sian beam (fi = 2.9 x 10~* sr), in addition to the 
map data. A flux calibration multiplier of 1.5±0.15 
was applied to both Bolocat and the image mo- 
saics to correct a Vl .O pipeline error (see All and 
IGinsburg et al.|[2"013l for a fuU discussion). 

The BGPS data pipeline removes atmospheric 
signal using a principle component analysis tech- 
nique that discards time-stream signals correlated 
spatially across the bolometer array. This effec- 
tively acts as an angular filter, attenuating angular 
scales comparable to or larger than the array field 
of view (see All, their Fig. 15). The implication 
is that the BGPS is not sensitive to scales larger 
than 6'. The effective angular size range of detected 
BGPS sources therefore corresponds to anything 
from molecular cloud cores up to e ntire clouds de- 
pend ing on the heliocentric distance ([Dunham et al.l 
1201 If) . In this work we refer to BGPS objects as 
"molecular cloud clumps" for simplicity, but recog- 
nize that distant sources are likely larger structures. 

2.2. Spectroscopic Follow-Up of BGPS Sources 

Several spectroscopic follow-up programs have 
been conducted to observe BGPS sources in a va- 
riety of molecular emission lines that trace the 
dense gas associated with molecular cloud clumps. 
These surveys provide both kinematic and chem- 
ical information, and are typically beam-matched 
to the BGPS to facilitate comparison to the dust- 
continuum data. From these observations, a line-of- 
sight velocity (i'lsr) "^^^ successfully fitted for each 
of more than 3,500 detected sources. A summary of 
spectroscopic programs is presented in Table [TJ 

* Available through IPAC at 
http : //irsa . ipac . caltech . edu/data/BOLOCAM_GPS 



In a pilot st udv (iSchlinginan et al.l 1201 ID and 
complete survey ([Shirlev et al.l2013D . all 6,194 Bolo- 
cat objects at £ > 7? 5 were observed using the 
Heinrich Hertz Submillimeter Telescope (HHT) on 
Mt. Graham, Arizona. These studies simultane- 
ously observed the J=3— 2 rotational transitions 
of HCO+ {v = 267.6 GHz) and N2H+ {ly = 
279.5 GHz). Because these molecular transitions 
trace fairly dense gas {nc-g « 10^ cm~'^j3, the line- 
of-sight confusio n seen in CO studies is largely ab- 
sent. In fact, iShirlev et al.l find only 2.5% of 
HCO''" detections have multiple velocity compo- 
nents. These objects, likely an overlap of two or 
more molecular cloud clumps along the line of sight, 
are not used in this study. Detectability in HCO"'" is 
a strong function of millimeter flux density, and the 
dete ction rate for the f ull HHT survey was « 50% 
fsee lShirlev et al.ll2013l for fuU details). Velocity fits 
to HGO"*" spectra constitute the bulk of the kine- 
matic data used in this study (N2H+ spectra were 
not used because the complex hyperfine structure of 
its transitions makes it difficult to fit ^lsr)- 

As a companion to the HHT observations, a 
subset of 555 BGPS sources were observed in the 
J=2-l rotational transition of CS {f = 97.98 GHz) 
using the Arizona Radio Observatory 12m telescope 
on Kitt Peak (Y. Shirley 2 012, private communica- 
tion; see iBallv et al.ll2010D . This subset was con- 
fined to 29° < £ < 31°, a region with a high density 
of sources looking toward the Molecular Ring and 
the end of the long Galactic bar. This transition of 
CS traces lower density gas (rioft « 5 x 10"^ cm^'^) 
than HC0+(3 - 2), and was detected in 45% of 
sources not detected by the HHT survey in this re- 
gion. 

Seeking to char acterize the physical properties 
of BGPS sources, iDunham et all (j201in used the 
Robert F. Byrd Green Bank Telescope to observe 
the lowest inversion transition lines of NH3 near 24 
GHz. They observed 631 BGPS sources in the inner 
Galaxy. The NH3 (1,1) inversion is the strongest 
ammonia transition at the cold temperatures of 
BGPS sources (T « 20 K), and we used this transi- 
tion exclusively for the NH3 velocity fits. 

2.3. The Spitzer GLIMPSE Survey 

Th e Spitzer GLIMPSE sur vey (jBeniamin et all 
I2003t IChurchwell et all I2OO90 was used to iden- 
tify mid-infrared extinction features associated with 
BGPS detected sources. The GLIMPSE survey area 
completely encompasses the BGPS for \b\ < 1?0 and 
i < 65° (there are several sections of the BGPS that 
flare out to |6| < 1?5, see All). We used the V3.5 

IRAC Band 4 mosaicfH {K = 7.9 ^m) to iden- 
tify absorption features. Point sources (stars) iden- 

'■* The effective density required to produce line emission 
with a brightness temperature of 1 K; may be up to several 
order s of magnitude smaller than the critical density l l Evans! 
[T999IV 

^^ Data product manual: 
http : //irsa . ipac . caltech . edu/data/SPITZER/GLIMPSE/ 
doc/glimpse l_dataprod_v2 . . pdf 
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TABLE 1 

Spectroscopic Follow-up Observations of Inner Galaxy 

BGPS Sources 



Species 


Transition 


V 

(GHz) 


ResoL^" 

{") 


(cm-3) 


N "^ 


Rcf. 


HCO+ 


J = 3-2 


267.6 


28 


10* 


6194 


1 


N2H+ 


J = 3-2 


279.5 


27 


10-1 


6194 


1 


CS 


J = 2-1 


97.98 


64 


5 X 10^ 


553 


2 


NH3 


(1,1) 


23.69 


31 


10^ 


631 


3 



References. — (1) IShirlev et all (IMl): (2 ) Y. Shirley 
(2012, private communication'): (3) IDunham et al.l l|2011f ) 

^ Beam FWHM 

^ Approxima te efTective density for line excitation at T = 20 K 

HEvanjI 19991 ') 

"^ Number of unique BGPS sources observed in this line 



tified in the Band 1 mosaics (Af, = 3.6 //m) were 
removed from the Band 4 images to accentuate dif- 
fuse emission (see ij4.2p . Stars were modeled as 
Gaussian peaks since the mosaicing process from 
individual IRAC frames produces a spatially vari- 
able PSF, hampering star-subtraction. The Band 4 
mosaics have an angular resolution ^ 2", and a 
pixel scale of 1'.'2. GLIMPSE images have under- 
gone zodiacal light subtraction based on a zodiacal 
emission model (see the data product manuaP^) . 
so signal remaining in the mosaics is Galactic in 
nature. There is, however, a significant effect due 
to scattering of light within the IRAC camera that 
causes the surface brightness of exten ded emission 
to ap pear brighter than it actually i^ (|Reach et al.1 
l2005f ). The method used in this study to correct for 
scattered light is described in ii4.21 and a deriva- 
tion of the correction factors required for quanti- 
ties measured from the publicly-available GLIMPSE 
mosaics is given in Appendix |B] 

3. DISTANCE PROBABILITY DENSITY 
FUNCTIONS 

3.1. Approach and Utility 

We introduce an automated distance determina- 
tion technique for molecular cloud clumps that al- 
lows for the joint application of many individual 
distance estimation methods. Bayes' Theorem pro- 
vides a framework for creating distance probabil- 
ity density functions (DPDFs) for dust-continuum- 
identified molecular cloud clumps that encode the 
confidence in source distances. Kincinatic distances 
derived from w^sr ^^d a Galactic rotation curve con- 
stitute the likelihood functions in the Bayesian con- 
text. Because these likelihoods are subject to the 
KDA, prior DPDFs based on ancillary data must 
be applied to constrain the distance estimates. The 
posterior DPDF is simply the product of the likeli- 
hood with the priors, suitably normalized. Relative 
amplitudes of the posterior DPDF at each distance 
along the line of sight (dg) correspond to the prob- 
ability of the source being at that distance. 

^^ See §4.11 of http://irsa.lpac.caltech.edu/ 
data/SPITZER/docs/lrac/lraclnstrumenthandbook/ 



Within this framework, any number of prior 
DPDFs may be applied to constrain the distances to 
molecular cloud clumps. This paper describes two 
such priors. The first, applicable to all molecular 
cloud clumps, is based on the Galactic distribution 
of molecular hydrogen. Because the scale height 
of the nrolecular disk is small, this prior favors the 
near kinematic distances for objects at high Galac- 
tic latitudes. The second prior involves the use of 
EMAFs. Not all molecular cloud clumps are vis- 
ible as absorption features, however, so this prior 
(described in detail in Q applies only to a subset 
of objects. To expand the collection of molecular 
cloud clumps with well-constrained DPDFs, addi- 
tional techniques {e.g. HISA, NIREX, etc.) would 
need to be applied. 

Not only do DPDFs provide a structure for apply- 
ing multiple techniques for distance discrimination, 
they also encode the distance uncertainty and level 
of confidence in the KDA resolution. When used 
to derive the mass or other property of a molecu- 
lar cloud clump, DPDFs provide a means for de- 
termining the associated uncertainty. The DPDFs 
derived in this work are computed out to a helio- 
centric distance of 20 kpc in 20-pc intervals. To fa- 
cilitate the use of integrated probabilities, DPDFs 
are normalized to unit total probability such that 
/,°°DPDFd(dJ = l. 

3.2. Extracting a Distance from the DPDF 

The proper use of DPDFs for calculating derived 
quantities is to build a distribution by randomly 
sampling distances from the DPDFs in a Monte 
Carlo fashion, preserving all information about dis- 
tance placement and uncertainty. There are appli- 
cations, however, that benefit from or require a sin- 
gle distance estimate with uncertainty (such as dis- 
tance comparisons with other studies). There are 
two primary distance estimates that may be derived 
from a DPDF. The maximum-likelihood distance 
(c^^l) is the distance which maximizes the DPDF. 
This represents the single best-guess at the distance 
for cases where a large fraction of the total proba- 
bility lies within a single peak. The associated un- 
certainty may be defined as the confidence region 
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around d^^ that encloses at least 68.3% of the in- 
tegrated DPDF, and whose limits occur at equal 
relative probability. This so-called isoprobability 
confidence region is generally asymmetric, and may 
represent lopsided error bars several kiloparsecs in 
size if both kinematic distance peaks are required 
to enclose sufficient probability. The full width of 
this uncertainty (FWgs), therefore, provides a direct 
measure of how well constrained a distance estimate 
is. Error bars produced in this way should not be 
considered Gaussian, as the 95.5% and 99.7% iso- 
probability confidence regions may be similar in size 
to the 68.3% error bars, or be radically different. 

An alternative single-value distance estimate is 
the weighted average distance (d), the first moment 
of the distribution, 



d„ DPDF d(d ) . 



(1) 



If the DPDF is well-constrained to a single peak, 
djj^j^ and d will be nearly equivalent. In cases where 
the KDA resolution is not well-constrained, how- 
ever, these distance estimates may be substantially 
different and d is not a good estimator of the dis- 
tance. The uncertainty associated with d may be 
computed from the second moment of the DPDF as 



/ d^2j3PDFd(dJ 
Jo 



d 



1/2 



(2) 



The a- represent the variance of the DPDF, and 
only approximate Gaussian confidence intervals for 
single-peaked DPDFs. Ultimately, the choice of a 
single-value distance estimate will depend on the 
specifics of the application; various cases are dis- 
cussed in ij6.1.2l 

3.3. Using DPDFs to Estimate Physical 
Parameters 

While distances to objects are often interesting 
in isolation, their primary use is to convert obser- 
vational quantities into physical properties of the 
object. DPDFs offer a simple way to propagate 
the uncertainties in distance through these calcula- 
tions. For example, the maximum-likelihood mass 
of a molecular cloud clump can be estimated as 



M,, 



S, , d. 



(3) 



where 5"^ ^ is the A = 1.1 mm flux density, and a 
contains the dust physics and temperature. Adop- 
tion of a DPDF representation allows marginaliza- 
tion over distance to obtain the expectation value 
of the mass: 



(M) = 



5,, d 2 DPDF d(dj 



(4) 



Practically, this integration can be accomplished by 
Monte Carlo methods, drawing a large number of 
distance samples from the DPDF and evaluating 
the average mass. Uncertainties in the expectation 
value can be determined using methods paralleling 
those used for distance above. 



Bimodal DPDFs again lead to complications, as 
the expectation value will commonly be found at 
a value with low probability. A maximum likeli- 
hood distance can be adopted to avoid this aesthetic 
feature, but marginalization over the distance re- 
mains the most rigorous approach. Ideally, addi- 
tional prior DPDFs should be applied in order to 
minimize bimodality. 

3.4. Kinematic Distance DPDFs 

Kinematic distances form the foundation for the 
Bayesian approach to distance estimation, com- 
puted from the intersection of the Galactic rotation 
curve projected along the line of sight, v^d^), with 
the observed molecular line Wlsr- Transformation 
of velocity uncertainties onto the distance axis is fa- 
cilitated by the use of two-dimensional probability 
density functions, -P(Wlsr; "^o)- 

The rotation curve function, i^rotc (^^lsr ; '^o)j is 
constructed as 



-ProtcKsR.c^o) =exp — 



v{d^)Y 



2ai 



(5) 



where the uncertainty cTvir is the magnitude of ex- 
pected virial motions within regions of massive- 
star formation, accounting for peculiar motions of 
individual molecular cloud clumps {— 7 km s~^; 
iReid et all |2009| PI. The function is Gaussian in 
Ulsr, and is centered along v{df^); if integrated over 
Ulsr, a uniform DPDF is obtained. The probabil- 
ity density function from spectral line information 
(^spcc) is a Gaussian centered at the measured wiine; 
with observed linewidth crfjjj,,, independent of dg. 
As with Protc, this function yields a uniform DPDF 
when integrated over w^sr . Since Protc does vary as 
a function of d^ , localized peaks in the (Wlsr , d^ ) 
plane result when it is multiplied by Pspoc- The de- 
sired one-dimensional DPDFkin is obtained by sub- 
sequent integration over w^^gj^ . 

DPDFkin is double-peaked and symmetric about 
the tangent distance for objects with i?gai < i?Oi 
and single-peaked otherwise. The iifd^) were com- 
puted us ing the flat rota tion curve of IReid et all 
i|2009D . ISchonrich eFaLl (|2010D subsequently de- 
rived newer estimates of the Solar peculiar motion, 
affectin g rotation c urve fits to the maser parallax 
data of IReid et al.l The updated values used here 
are Rq = 8.51 kpc, and Oq = 244 km s"^ (M. Reid 
2011, private communication). The new solar mo- 
tion values also had the effect of decreasing the mag- 
nitude of the apparent Galactic counter-rotation of 
high-mass star forming regions, an effect likely aris- 
ing from molecular gas interacting with the spiral 
potential, from 15 km s~^ to 6 km s^^. 

Kinematic distances are sensitive to the slope of 
w(dQ), itself a function of Galactic longitude. For 

^^ This is the expected virial velocity, per coordinate, for 
an individual object (i.e. molecular cloud clump) within a 
high-mass s tar-forming regio n of mass ~3xl0* Mq and ra- 
dius ~1 pc l|Reid et al.i[2ro9l ). 
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lines of sight along & w 0° within ^ 10° of the 
Galactic longitude cardinal directions, w(rfQ) is ei- 
ther very flat or sharply peaked; small departures 
from circular motion therefore translate into large 
deviations in derived kinematic distances. Further- 
more, since v^d^) is derived assuming circular or- 
bits about the Galactic center, radial streaming mo- 
tions of the gas are not accounted for, meaning 
that DPDF-derived distance estimates carry the ba- 
sic limitations of any kinematic distance determina- 
tion. To minimize the effects of non-circular motion, 
regions known to have significant streaming must 
be excluded from consideration. In particular, the 
prese nce of the long Galactic bar at _Rgai 5. 3 kpc 
(|Fuxl 119991: iRodriguez-Fernandez fc Combed [2008[) 
and its associated radial streaming motions restrict 
the use of kinematic distance measurements to loca- 
tions outside this radius. In the Galactic longitude- 
velocity {i — v) diagram, these restrictions amount 
to excluding much of \i\ < 20°. Features at low lon- 
gitude known to be outside the Galactic bar (such 
as the Scutum-Ce ntarus arm, also labeled as the 
"Molecular Ring": lDame et al.|[200l their Fig. 3), 
may be considered to have roughly circular orbits, 
and are included in this study. 

3.5. Prior DPDFs for Kinematic Distance 
Discrimination 

Prior DPDFs are required to discriminate be- 
tween the kinematic probability peaks for objects 
within the solar circle. DPDFkin is symmetric about 
the tangent point, so prior DPDFs based on an- 
cillary Galactic plane data must be asymmetric to 
provide useful distance constraints. 

The Galactic distribution of molecular gas serves 
as an envelope inside which molecular cloud clumps 
may form. The prior DPDFhs is defined to be pro- 
portional to the volum e density from the m olecu- 
lar hydrogen model of iWolfire et al.l (|2003D along 
a line of sight. This model consists of a Molec- 
ular Ring component with a decaying exponential 
toward the outer Galaxy; the vertical distribution 
is Gaus sian with a half-width at half maximum of 
60 pc (jBronfman et al.l 119880 , flaring outside the 
Solar Circle. While this distribution is symmetric 
about dtan along the Galactic midplane, the nar- 
row vertical extent of the molecular layer sets a 
strong prior on higher-latitude objects. The relative 
amount of H2 beyond the tangent point for lines of 
sight at l&l > 0?3 is small, generating the needed 
asymmetric function for molecular cloud clumps at 
larger Galactic latitude. 

The prior DPDF based on EMAFs was computed 
from a pixel-by-pixel morphological matching be- 
tween millimeter dust-continuum emission and mid- 
infrared dust absorption features. The derivation of 
DPDFomaf is described in detail in the next section. 

4. INFRARED-MILLIMETER 
MORPHOLOGICAL MATCHING 

Morphological matching is based on the com- 
parison between synthetic 8-/xm images computed 



from millimeter flux density measurements and 
GLIMPSE 8-/im maps processed to match the an- 
gular resolution of the BGPS. This section describes 
the creation of both the synthetic and processed 8- 
//m images, as well as the mechanics of computing 

DPDFcnaf. 

4.1. Creation of Synthetic S-jim Images 
4.1.1. Radiative Transfer Assumptions 

Creation of synthetic 8-/xm images explicitly as- 
sumes that the dust seen in emission in the BGPS 
is the same dust that extincts mid-infrared light. 
When converted into a mid-infrared optical depth, 
BGPS observations represent dark clouds which 
may be placed at different heliocentric distances 
within a model of diffuse Galactic 8-/im emission. 
A series of synthetic images generated in this man- 
ner were compared with mid-infrared observations 
to compute the DPDFgmaf- 

We assumed a simple radiative transfer model 
to describe the observed mid-infrared intensity ab- 
sorbed by a cold molecular cloud clump immersed in 
a sea of diffuse emission (assuming that the absorb- 
ing cloud has no emission). The intensity observed 
within an EMAF (/cmaf) is 



^cmaf ^^ -'back ^ 



/) 



fore 



(6) 



where /back and /fore are the background (from the 
cloud to large heliocentric distance) and foreground 
(between the observer and the cloud) intensities, re- 
spectively, and Tg is the mid-infrared optical depth 
of the cloud. The total intensity along a line-of-sight 
in the absence of absorption is I^^^^ = /back + /fore- 
Defining the fraction of the total intensity that lies 
in front of the cloud as /fore — Iiotc/Imiti allows 
Equation ^ to be written as 

/cmaf = [(1 - /fore) e^^'s + /fore] /mir • (7) 

This parameterization frames the observed EMAF 
intensity in terms of decrements below the un- 
extincted intensity in the vicinity, and provides the 
basis for creating synthetic 8-/im images. It fol- 
lows quickly from Equation ([7]) that clouds optically 
thick in the mid- infrared (r^ ~ 1) will still have a 
10% difference between /emaf and /^^j^ {i.e. easily 
detectable) for /fore as large as 0.85. Calculation of 
Tg and /fore are described below, and the estimation 
of 4„j, from GLIMPSE data is discussed in gjl 

4.1.2. 8-fini Optical Depth from the Millimeter Flux 
Density 

The mid-infrared optical depth of an EMAF can- 
not be measured directly from the GLIMPSE mo- 
saics without significant assumptions, but it may be 
estimated from millimeter data. Thermal dust emis- 
sion is optically thin at millimeter wavelengths, so 
the observed BGPS fiux density {S-^ ^ ) may be writ- 
ten as 

Si.i = B^^{Td) Tj 1 rigQps , (8) 

where B^ i(/d) is the Planck function evaluated at 
A = 1.1 mm and dust temperature T^, and J^bgps — 
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2.9 X 10~^ sr is the solid angle of the BGPS beam. 
The millimeter optical depth (r^ ^ ) was computed 
assuming the dust opacity (k^ ^) for grains with thin 
ice mantles, coagulatin g at 1 0^ cm~'^ for 10^ years 
([Ossenkopf & Henning 1994 Table 1, Column 5; 
called 0H5 dust). Interpolation of 0H5 dust opaci- 
ties to the central frequency of the BGPS bandpass 
yields k^ ^ = 1.14 cm^ g~^ of dust (All). A molecu- 
lar cloud clump with r^ j^ = 10~^, which corresponds 
to (S'j J ~ 0.9 Jy), has a beam-averaged molecular 
hydrogen column density Ki 2 X 10^^ cm~^. 

The 8-/im optical depth is related to t^ ^ by the 
ratio of the dust opacities in the two bandpasses, 
i?K = Ks/i^i 1 • We calculated the mid-infrared dust 
opacity by assuming a d ust emission spectr um in- 
cluding PAH molecules (jDraine fc Lill2007[) . find- 
ing the average attenuated intensity across IRAC 
Band 4, and extracting a band-averaged opacity 
Kg = 825 cm^ g~^ of dust (see Appendix IX|). At 
the 33" resolution of the BGPS, the beam-averaged 
8-/im optical depth is therefore 
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Rk 



13.0K/Td 

: 0.778 ' 



I \ f s, 



gl3. OK/20. OK _ I J I 1 Jy 



(9) 



The function T{Td) has units of inverse flux density, 
and is normalized to 20 K in Equation (|9]). Because 
Tg is a function of Rk [i-e. both millimeter- wave 
emission and mid-infrared absorption depend only 
on the dust), the dust-to-gas ratio is not relevant 
to the distance estimation method. Owing to the 
nearly three orders of magnitude difference in dust 
opacity between the millimeter and mid-infrared, a 
value of Tg = 0.1 corresponds to a column of only 
JV(H2 ) pa 3 X 10^^ '^'^.~^ assuming Afg^jMv » 
0.05 (jlndebetouw et al.|[2005t iRoman-Zuiiiga et al.l 
I2007D . Therefore, molecular cloud clumps with col- 
umn densities > 10^^ cm~^ will be mostly opaque 
at A = 8 /xm. 

Using Equation ([9]) to obtain an H-fim op- 
tical depth requires a dust temperature (T^). 
Since we are ignorant of Td within each molecu- 
lar cloud clump used in this study, we assumed 
that all source s are at the same temperature. 
iBattersbv et al.l (|2011[ ) showed that mid-infrared- 
dark molecular cloud clumps generally span the 
temperature range 15 K < T^; < 25 K. Therefore, 
Td = 20 K is a reasonable representation for BGPS 
sources as a group. Variation of the assumed Td 
affects the KDA resolutions for some sources, and 
is discussed briefly in i JG.l.ll With molecular cloud 
clump dust temperatures derived from Herschel Hi- 
GAL data, more precise DPDFs for individual ob- 
jects may be derived using the present methodology. 

4.1.3. S-fim Foreground Fraction from a Galactic 
Emission Model 

Absorption features seen at A = 8 fira are as- 
sumed to be the result of dense clouds immersed 
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Fig. 1. — Galac tic mid-infrared emissi on model computed 
with Hyperion l IRobitaille et al.l 120121 ) viewed from the 
North Galactic Pole. The model, viewed through the IRAC 
8.0-/^m bandpass, is shown on an inverted square-root in- 
tensity scale. The Sun is located at {x,y) = (0,-8.5 kpc), 
and solid diagonal lines represent the limits of the GLIMPSE 
survey {\£\ = 65°). The dashed line marks the low-latitude 
{£ = 7? 5) limit of this study. DPDFs were computed out to 
rf„ = 20 kpc (curved contour). 



in a smooth emission distribution, punctuated by 
regions undergoing active star formation. While 
small-scale structures are difficult to model, the 
broader diffuse emission is a more tractable prob- 
lem. Creation of synthetic 8-/im images via Equa- 
tion ([7]) requires a three-dimensional model for the 
Galactic 8-fim emission distribution. 

The recent numer ical Galactic ste ll ar an d dust 
emission model of iRobitaille et al.l (|2012l here- 
after R12), computed using the Monte-Carlo three- 
dimensional rad iative transfer code HyperioisF^ 
(jRobitailld 1201 If ), offers a self-consistent estimate 
of diffuse Galactic emission that is well-matched to 
observed quantities. We used the final model pre- 
sented in R12, whose parameters were chosen to fit 
the Galactic latitude and longitude intensity dis- 
tributions from seven bandpasses in the mid- to 
far-infrared. This model features two major and 
two minor spiral arms with Gaussian radial pro- 
files, a lack of dust in the inner few kiloparsecs of 
the Galactic disk (dust hole; correlated with the 
dearth of molecular gas in this region) , and a mod- 
ified PAH abundance rel ative to the favored model 
from iDraine fc Lil (j2007| ) . An analysis of the con- 
tributions from various stellar populations and dust 
grain sizes to the total intensity in each bandpass 
indicates that some 96% of the emission detected in 
IRAC Band 4 images comes from PAH molecules 
(R12). 

A three-dimensional (£, b^d^ ) data cube of 8-^m 

http://www.hyperion-rt.org 
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TABLE 2 
Comparison of Hyperion Model Parameters 



Category 


Parameter 


R12 




This Work 


Grid^ 


Nr 

N^ 
N^ 

^Imax (pc) 


200 

100 

50 

3000 




200 

200 

44 

1000 


Wavelength'' 


N bins 
Range {fJ.m) 


160 
3 < A < 140 


22 
6 < A < 10 


Image'' 


Observer i?gai (kpc) 
Observer z (pc) 
Longitude Range (°) 


8.5 
+15 
65 > £ > ^ 


-65 


8.5 

+25 

65 > £ > -65 



^ Af = number of grid cells in this dimension 
^ Wavelengths at which the model images were computed later 
convolved with instrument bandpasses to create simulated obser- 
vations. 
•^ Parameters related to observer within the grid. 



emission was generated from the radiative transfer 
code using model grid and ima ge parameters sli ghtly 
modified from those used by iRobitaille et all Ta- 
ble [2] lists the comparison of Hyperion input pa- 
rameters between R12 and the present study. Pri- 
mary differences include an increase in azimuthal 
resolution of the cylindrical grid, a restriction of 
the vertical extent of the grid to \z\ < 1 kpc (to 
match the region of the Galactic plane probed by 
the latitude range |6| < 1°), and limiting the wave- 
length range used in computing output images. The 
model was computed within a box 30 kpc on a side, 
containing the entire modeled stellar disk (R12); a 
face-on view of the model Milky Way as seen from 
the north Galactic pole is shown in Figure [1] Lines 
of sight out to 20 kpc (the distance used for DPDF 
generation) lie entirely within the simulation box 
for 1^1 < 48°. Beyond this longitude, however, the 
edge of the box retreats to only d^ « 16.5 kpc by 

\e\ = 65°. 

A series of {£, b) images of the Galactic plane 
containing only emission between the observer and 
some distance di were computed using a useful Hy- 
perion feature, using a resolution of 3' in latitude, 
and 15' in longitude for computational reasons. Im- 
ages were computed for each of 22 wavelength bins 
logarithmically spaced from 6 to 10 fim (closely 
matching the wavelength bins used by R12 for this 
part of the spectrum), and were then convolved 
with the IRAC Band 4 transmissio n curve to yield a 
single simulated Spitzer image (see IRobitaille et al.l 
120071 for complete details) . The three-dimensional 
image was constructed by stepping di outward in 
100-pc intervals and depicts the cumulative 8-fiin 
emission out to each di. The resulting cube com- 
prises 200 steps, with the final image slice including 
all model emission to the edge of the box, equivalent 
to the collapsed profiles presented in R12. 

Hyperion cannot treat sources individually, but 
rather uses "diffuse" sources of emission in each 
grid cell. These diffuse sources are generated from 
the probability of emission from various popula- 
tions of stars and similar objects {e.g. planetary 
nebulae, H II regions), assigned a spectrum cor- 



responding to the appropriate spectral class, and 
given a total luminosity based on the number of 
"real" sources the cell represents. While most of 
the emitting populations have a smooth spatial dis- 
tribution, relatively rare sources with concentrated 
emission at A = 8 /im (such as H II regions) are 
sprinkled throughout the box according to the un- 
derlying stellar distribution model. Very nearby ob- 
jects (dg < 0.5 kpc) appear quite bright, and cause 
"hot-pixel" effects in the computed images of the 
Galactic plane. These objects blend into the back- 
ground for images computed from large Galactocen- 
tric position (e.g. Fig. [Ij, or are averaged out in 
collapsed longitude or latitude distributions (R12). 
To ameliorate the effect of these objects in the com- 
puted {i, b) images, we ran seven realizations of the 
model, each with a different random-number seed, 
then median-combined the realizations of each di 
slice. Since the underlying distribution of sources 
is fixed, nearby bright sources often appear in the 
same pixel in the output images; the number of re- 
alizations was chosen to be large enough such that 
median combining the realizations removes most of 
these outliers. To eliminate any remaining outliers 
and reduce noise, the combined (£, b) images were 
median smoothed with a 3 pixel x 3 pixel box. 

The foreground fraction was computed from the 
intensity cubes by dividing each (£, b) image slice by 
the final slice. The final fits data cubes of 8-/im 
intensity and /foro for both the northern and south- 
ern Galactic plane (|£| < 65°) are publicly available 
with the BGPS archive. To illustrate the Galac- 
tic features present in the modeled cube, ffore{(^,d^) 
for the Northern plane along 6 = 0° is shown in 
Figure [21 with contours and grayscale representing 
its value from to 1. Since PAH molecules con- 
tribute the bulk of the model emission, the dust 
hole towards low longitude is visible as a flattening 
of /fore(rf0)- The Molecular Ring / Scutum tangent 
at ^ w 30° appears where /fore grows quickly as 
a function of distance. The limited distance range 
caused by the model box size is represented by the 
1.0 contour for i > 48°. The tangent distance as 
a function of longitude (black dashed line) spans 
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Fig. 2. — Foreground fraction of Galactic 8-fiai emission 
in the northern Galactic plane derived from the Hyperion 
model as a function of (£,d„) along 6 = 0°. Grayscale 
and contours represent /forei with the unlabeled 1.0 contour 
marking the edge of the box in Fig. [T] for £ > 48°. The thick 
black dashed line follows the tangent distance as a function 
of Galactic longitude, and the vertical dot-dashed line marks 
the £ = 7? 5 lower limit of this study. 

the range 0.45 < /fore ^ 0.6, implying that clouds 
that are optically thick in the mid-infrared should 
be visible beyond dtan- 

4.1.4. Computing the Synthetic Images 

Synthetic images (/cmaf) for a given BGPS ob- 
ject are computed using Equation ^. The optical 
depth is modeled as a two-dimensional image, con- 
structed by applying Equation ^ to the BGPS map 
data. The estimate of the total mid-infrared emis- 
sion (/mjr) is also a two-dimensional image, and its 
creation is discussed below. Because of the coarse 
resolution of the /fore model, we simply extracted 
the one-dimensional fioreid^) at the (€,6) of the 
BGPS object. The combination of these elements 
yields a cube of synthetic data to be compared with 
the processed GLIMPSE images. 

4.2. Processing of GLIMPSE 8-^m Images 

Mid-infrared properties of dust-continuum- 
identified molecular cloud clumps were derived 
from the Spiteer/GLIMPSE mosaics. Further pro- 
cessing of these images was required to estimate the 
total mid- infrared intensity (/mir) i'^ the vicinity 
of an EMAF, and to produce a smoothed, star- 
subtracted map, containing features and angular 
scales comparable to (sub-)millimeter data. The 
example source G035.524-00.274 (BGPS #5647) is 
used to illustrate the processing products in Fig- 
ure [21 The first step was to remove individual stars 
because they contaminate estimates of broader 
diffuse emission and (sub-)millimeter observations 
are not sensitive to them. Star locations were 
identified by searching for bright, unresolved 

objects in th e 3.6-Atui mosai cs using DAOFIND 

(jStetsonl 119871: ILandsmanl |1995l ) with a threshold 
of 20 MJy sr~^. A Gaussian was fit to the 8-/im 



image at the location of each identified star, then 
subtracted. This method of star subtraction was 
deemed optimal because PSF variations across the 
survey mosaics meant that PSF-based approaches 
could not be applied. Star subtraction in this 
manner did, however, leave clear low-level residuals 
(Fig. ^p). Since later processing smooths the 
resulting images to the BGPS resolution, residuals 
are largely unimportant. However, to ensure that 
poor star subtraction or other effects did not effect 
distance estimation, by-eye evaluation of each 
potential EMAF for contamination was performed. 

For further processing of the GLIMPSE data, 
6' X 6' postage-stamp images were extracted from 
the star-subtracted 8-/xm mosaics for each Bolocat 
source. These postage stamps, centered on the lo- 
cation of peak millimeter flux density, limit consid- 
eration of mid-infrared variations to the immedi- 
ate vicinity of a molecular cloud clump in addition 
to providing computational expediency. The first 
postage-stamp image created for a given BGPS ob- 
ject is a version of the star-subtracted GLIMPSE 
mosaic, re-pixelated and aligned to the 7'.' 2 scale 
of the BGPS images. This image was used to en- 
sure that locally bright emission did not interfere 
with derived mid-infrared intensities, and to de- 
termine the likely intensity range containing I^^^ 
around the object. A pixel intensity histogram of 
the image was constructed with 1 MJy sr~^-wide 
bins, and the background was defined as intensi- 
ties within its full-width at half maximum. Pixels 
within an 8' x 8' section of the native- resolution star- 
subtracted GLIMPSE mosaic having intensities in 
the defined range were used to fit a quadratic sur- 
face using a linear, least-squares optimization. This 
surface, repixelated and scaled as above, comprises 
the postage-stamp estimate of /„jp. (Fig. [31i). This 
estimate of background pixels ignored high pixel val- 
ues from star residuals and low pixel values from 
EMAFs. 

The IRAC camera on Spitzer suffers from inter- 
nal scatt ering of light which affects instrument cal- 
ibration ([Reach et al.l 120051 ) . Point-source photom- 
etry is unaffected by the scattering due to the cali- 
bration technique employed, but extended emission 
(such as the Galactic plane) will appear brighter 
due to scattering into each pixel. Gorrecting for 
this effect should be done on a frame-by-frame ba- 
sis, but was not accounted for in the GLIMPSE 
pipeline (S. Carey 2010, private communication). 
To approximately correct for the scattering, an es- 
timate of the scattered light was subtracted from 
the postage-stamp images for each BGPS source. 
The postage-stamp size was chosen to be near the 
5f2 X b'.2 FOV of IRAC, and the I^^^ fit serves 
as the estimate of the light available to be scat- 
tered within a single IRAC frame. This estimate is 
only approximate, as the /^jj^ fit explicitly excludes 
very bright and very dim emission within a frame; 
for frames with regions of bright emission, the de- 
rived correction factor will be a lower limit, and vice 
verse for frames containing extensive dark clouds. 
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Fig. 3. — BGPS and processed GLIMPSE data for example object G035. 524-00. 274. All panels are 6' X 6' postage-stamp 
images (see text), and the pink circle identifies the 40" top-hat BGPS equivalent aperture, centered on the location of peak flux 
density, (a) Cutout of the star-subtracted GLIMPSE image at native resolution. Cyan ellipses mark IRDCs identified in the 
IPeretto &: Fulled l|2009l ) catalog; note that the dark cloud associated with BGPS source G035. 478-00. 298 (lower-right corner) is 
not included in that catalog. (6) BGPS map data with Bolocat source boundary (light cyan), (c) Star-subtracted GLIMPSE 
cutout smoothed to 33" and resampled to 7" 2 pixels to match the BGPS maps. Color contours represent logarithmic flux 
density levels from BGPS (Jy beam"^). (d) Estimate of the total mid-infrared intensity, If^i^, as a quadratic surface fitted 
to background pixels as described in the text. Contours are drawn to show the variation in /^jjj^ over the postage stamp 
(MJy sr~^; ticks point to higher values), and the cyan aperture marks the region used to estimate {Ifjn^^) for this source (see 
i|5.1|l . The grayscale colorbar represents the common (logarithmic) intensity scale for all three GLIMPSE panels. (A color 
version of this figure is available in the online journal.) 



The infinite-aper ture intensity corre ction for IRAC 
Band 4 is 0.737 (jReach et al.ll2005[) . meaning that 
^ = 0.263 is the scattered light fraction. Assuming 
that /j^^jp, represents the total incident light, we sub- 
tracted ^ X median(/^jj^) from each postage-stamp 
image to remove scattered light. 

Reduction of the GLIMPSE angular resolution 
was necessary for direct comparison with the syn- 
thetic images created using Equation ^. Since 
bright emission in the vicinity is scattered into 
an EMAF, removal of the scattered light must be 
done prior to smoothing. The scattering-corrected 
extracted postage stamps were smoothed with a 



FWHM = 33" Gaussian kernel, then re-pixelated 
and aligned to match the BGPS images (Fig.[3]c). 

4.3. Morphological Matching 

Derivation of DPDFcmaf relies upon the com- 
parison of the (sub-)millinieter emission and mid- 
infrared absorption of dust in cold molecular cloud 
clumps. The series of synthetic images were 
matched against the smoothed GLIMPSE postage 
stamp images (Fig. [3]c) . For small d^ , the synthetic 
cloud appears darkest, without foreground light fill- 
ing in the absorption feature. At larger heliocentric 
distance, /foro grows, and the synthetic image con- 
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Fig. 4. — Morphological matching for example object G035. 524-00. 274. The millimeter source to the lower-right of the marked 
contour is a separate Bolocat object, (a) Synthetic 8-^m image calculated via Equation (O. (6) BGPS postage-stamp image, 
showing the restricted region used for the morphological matching (see text), (c) Smoothed GLIMPSE map against which 
the synthetic images are compared, (d) Prior DPDFs from the morphological matching (black) and molecular gas distribution 
(blue dot-dashed), and the posterior DPDF (red), including the kinematic distance likelihood (see text). The gray dotted line 
represents /fore extracted from the model cube along {i,b), and the green dashed line marks the tangent distance. (A color 
version of this figure is available in the online journal.) 



verges upon /^„^ (Fig. [31fi). 

The Galactic 8-/ini emission model of ^4.1.31 
describes smooth, diffuse emission against which 
EMAFs are visible, but actual Galactic emission is 
more complex. To match more closely the assump- 
tion of the model, the angular region over which the 
synthetic and observed sky are compared must be 
restricted. As the observational definition of a sin- 
gle molecular cloud clump, we began with a source's 
Bolocat contour delineati ng the maximu m extent 
of the comparison region (jRosolowskv et a l. 2010). 
Since the synthetic image can never be brighter than 
/„jj,, the matching process is adversely affected by 
bright mid-infrared emission in the vicinity of a 
BGPS source. To ameliorate this effect, pixels in the 
smoothed GLIMPSE postage stamp were excluded 



from the matching region if their value exceeded the 
corresponding value in the /„jj^ image. 

An overview of the morphological matching pro- 
cess is presented in Figure |4] for the same ob- 
ject as in Figure |31 The synthetic 8-^m image is 
shown in panel (a) for the distance which maxi- 
mizes DPDFemaf (sce bclow). Panels (6) and (c) 
are identical to Figure |3l except that the source 
contour now marks the restricted matching region 
due to bright mid-infrared emission on the perime- 
ter of the EMAF. Panels (a) and (c) are shown on 
a common linear grayscale to illustrate the match 
between the observed extinction and that predicted 
from thermal dust emission. The various DPDFs 
for this source are shown in panel (d), and are de- 
scribed below. 
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Quantification of the match as a function of dis- 
tance was accomplished by constructing a ^ statis- 
tic from a pixel-by-pixel comparison within the 
matching region. The estimate of the error in each 
pixel was derived from Equation ^ by propagating 
the uncertainty in the optical depth map as 



a.y„(^,6) = 4„,,(^,6)e--«(^^'')Ta3 



(10) 



where Equation (J9l) defines r^ and T, and a^ is 
the median absolute deviation of the BGPS postage- 
stamp image. This estimate of the uncertainty 
places more weight on the portions of the image 
with larger BGPS flux density. The statistic was 
computed for synthetic images at 100-pc intervals 
along the line of sight, yielding ^{d^^. 

A preliminary DPDFomaf was computed using the 
formal probability of the Ax^ statistic. The num- 
ber of degrees of freedom was taken as the inte- 
ger number of BGPS beams in the matching re- 
gion (A'pixois / 23.8 pixels beam~^; All) minus one, 
since only beam-scale structures are independent 
and distance is a fitted parameter. For most sources, 
the DPDFomaf bas a broad peak (several kilopar- 
secs wide), and falls sharply where the A(Xred) s^" 
ceeds unity. Because of the sharp cutoffs, it tends 
to very strongly favor one kinematic distance peak 
over the other. If the Galaxy truly consisted of dark 
molecular cloud clumps embedded within broad 
diffuse mid-infrared emission, this formulation of 
DPDFomaf would be appropriate. However, the 
Galaxy is punctuated with regions of stronger 8-/xm 
emission that violate the simple radiative transfer 
of Equation ©, and the DPDFomaf should contain 
a systematic uncertainty that allows non- vanishing 
probability at the non-favored kinematic distance 
peak. 

Experimentation with alternative approaches 
that allow a systematic uncertainty led to the se- 
lection of DPDFomaf fx (x^)~^, where /3 is a posi- 
tive scalar of order unity. This class of DPDFomaf 
have FWHM comparable to the formal probability, 
but greater width at low likelihood, and hence rarely 
goes to zero probability until far from the peak. The 
parameter /3 may be used to tunc the width of the 
function, with larger values leading to narrower dis- 
tributions. Since the sharp cutoff of the DPDFemaf , 
not the width of the peak, is what appears prob- 
lematic in light of complex Galactic emission, we 
selected /3 = 2 to reproduce the widths of the for- 
mal probability DPDF. To verify the validity of this 
choice, we computed the GRS distance matching 
success rate (see i j5.3.2p as a function of /3, and 
found no dependence on the width of DPDFomaf- 

The resuhing DPDFs for object G035. 524-00.274 
(BGPS #5647) are shown in Figure W-. The prior 
DPDFh2 (blue dot-dashed) favors the near kine- 
matic distance, since this line of sight looks out the 
bottom of the molecular layer. The gray dotted 
line shows the /foro('^Q) from the numerical model. 
The morphological matching process, represented 
by DPDFomaf (black solid) , could not make the syn- 



thetic image dark enough to match the smoothed 
GLIMPSE image, forcing the prior to peak at 
dg = kpc. The posterior DPDF (red) clearly 
reflects the distance selection, with the near peak 
containing > 95% of the integrated probability, al- 
though there remains some probability contained in 
the far kinematic distance peak at dg w 11 kpc. 

5. RESULTS 

5.1. EM AF- Selected Molecular Cloud Clumps 

5.1.1. Spatial and Kinematic Selection Criteria 

We derived posterior DPDFs for the subset of 
BGPS sources that have a measured v-^^^ from 
molecular spectroscopy, and are selected by the 
presence of an EMAF. Spatially, this set is defined 
by the GLIMPSE-BGPS overlap, limiting the upper 
end of the Galactic plane at ^ = 65?25 and a lati- 
tude spread of |6| < 1?0. Kinematic considerations 
restrict the regions of the £ — v diagram (Fig. [5]) that 
may be considered at lower Galactic longitude down 
to the i = 7°5 limit of the spectroscopic surveys. 

The colored image in Figure [S] is the latitude- 
integrated "'^^CO(l-O) intensity from l^ame et al] 
()200lD , and is shown as an indicator of molecular gas 
location and kinematic conditions. The presence of 
a lon g Galactic bar (i?bar '^ 4 kpc; iBeniamin et al.l 
|2005[) implies significant non-circular motion at £ < 
30°. Regions at these longitudes in the £—v diagram 
associated with the Molecular Ring feature (cf. 
Dame et "all 120011: iRodriguez-Fernandez fc Combed 



2008() . however, are likely at i?gai > 4 kpc. To in- 



clude the Ring but exclude bar-related gas, we disal- 
lowed the two hashed regions in Figure 5. The upper 
region is bounded by v^^^^^ = (3.33 km s~^)x^(°) -|- 
15 km s~^, and includes the higher- velocity gas in- 
side the Ring. The lower region excludes the 3- 
kpc expanding arm, and is bounded by v-^^-^ = 
(2.22 km s-^)xi{°) - 16.7 km s^^. Both regions 
are defined only for i < 21°. The upper hashed 
region does not extend past this point because the 
Molecular Ring feature extends to the tangent ve- 
locity at larger longitudes; the lower region is lim- 
ited because the 3-kpc ar m has its tangency here 
(jDame fc ThaddeusI l2008l ). There is likely over- 
lap between Ring objects with nearly circular mo- 
tions and objects in bar-related streaming orbits at 
21° < i < 30°, so kinematic distance estimates in 
this range, including those derived here, should be 
used with caution. 

Black circles mark the locations of the BGPS 
molecular cloud clumps for which a DPDFomaf was 
computed, and the histogram summarizes their lon- 
gitude distribution. Stars mark the masers used for 
distance comparison (see ij5.3.ip . 

5.1.2. Mid-Infrared Selection Criteria 

Automated classification of dust-continuum- 
identified molecular cloud clumps as EMAFs was 
achieved using the mid-infrared contrast computed 
from the smoothed GLIMPSE images at the loca- 
tion of the BGPS source. The peak contrast was 
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defined as 



C = \- 



(4„h) ' 



(11) 



where the intensity values were measured from the 
processed postage-stamp images described in §4.21 
for each Bolocat object. Due to the varied sizes 
and shapes of EMAFs, standardized intensities were 
measured in a 40" aperture around the location of 
peak BGPS flux density (pink circles in Fig. [3]). 
The value of /,„i,i is the minimum intensity within 
the aperture measured from the smoothed star- 
subtracted GLIMPSE image (Fig. E):;), and (4„^.) 
is the mean of the I^-^^ postage-stamp image within 
2' of the peak of millimeter flux density (cyan cir- 
cle in Fig. [Sji). A preliminary contrast threshold 
of C > 0.01 was implemented in the automated 
source selection to minimize the number of spuri- 
ous matches caused by unrelated variation in the 
GLIMPSE 8-/Ltm mosaics. This threshold also re- 
jects BGPS sources that are mid-infrared bright, as 
those objects have negative contrast. 

All molecular cloud clumps meeting the above 
selection criteria were examined by eye to ensure 
their suitability for deriving a DPDFomaf- Bolo- 
cat objects were not assigned a DPDFcmaf for the 
following types of deficiencies: (1) there was evi- 
dence of poor star subtraction contaminating /min; 
(2) there was very bright mid-infrared emission 



(^^un > 200 MJy sr"^) within 2' of the location of 
peak BGPS flux density that could bleed into the 
40" aperture or significantly affect the IRAC scat- 
tering correction; (3) the postage-stamp estimate 
of I^jpt was contaminated by excessive bright emis- 
sion or dark extinction; or (4) the morphology of 
dark regions in the GLIMPSE image clearly did not 
correspond to that of the millimeter emission. By- 
eye exclusion removed approximately 40% of sources 
meeting the initial automated selection criteria. 

Properties of rejected sources were analyzed to re- 
veal that nearly all very-low contrast sources were 
spurious matches (deficiency type 4, see above). Ad- 
ditionally, BGPS objects located in fields of locally 
very bright mid-infrared emission were almost all 
excluded from the final source list (types 2 and 
3). As a result, the contrast cutoff was increased 
to C > 0.05, and two additional automated selec- 
tion criteria were introduced. First, the restriction 
(-^MiR ) ^ 100 MJy sr~^ was placed to remove sources 
whose background estimate indicates significant dis- 
agreement with the 8-/xm-emission model, as large 
discrepancies may lead to improper distance esti- 
mates (type 3). Second, to automatically reject 
sources near very bright emission, the re-pixelated 
unsmoothed postage-stamp images were checked for 
pixels with /s^m > 200 MJy sr^^ within 2' of the 
image center; sources with more than 10 such (7'.'2) 
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TABLE 3 
Observed & Derived Properties of EMAF-selected BGPS Molecular Cloud Clumps 



BOPS VI. Catalog Properties'' 



Velocity Data 



Catalog 


I 


b 


540" 


''lsr 


Ref. 


Mid-Infrared 


KDA<= 


P <* 


•^ml'^ 


df 


Number 


n 


n 


(Jy) 


(kms-1) 




Contrast 


Resol. 




(kpc) 


(kpc) 


4638 


30.990 


0.329 


0.186(0.048) 


79.1 


2 


0.21(0.03) 


N 


0.88 


4.60tg1^ 






4639 


30.990 


0.385 


0.108(0.054) 


78.7 


2 


0.08(0.03) 


F 


0.91 


9.S8t°o.ll 






4650 


31.016 


-0.001 


0.280(0.057) 


74.5 


1 


0.25(0.04) 


N 


0.88 








4653 


31.026 


-0.113 


0.289(0.067) 


76.8 


1 


0.40(0.03) 


N 


0.83 


4 50+°-^^ 






4655 


31.032 


0.783 


0.267(0.110) 


51.0 


1 


0.38(0.05) 


N 


0.99 


3.24tg:^3^ 






4715 


31.226 


0.023 


0.381(0.077) 


74.5 


1 


0.40(0.02) 


N 


0.86 


4.4ot2:S 






4749 


31.342 


-0.149 


0.111(0.048) 


42.0 


1,3 


0.15(0.03) 


N 


0.91 


2.84tg-il 






4769 


31.432 


0.167 


0.117(0.039) 


101.7 


3 


0.06(0.05) 


U 


0.51 








4770 


31.436 


-0.103 


0.122(0.039) 


89.4 


1,3 


0.18(0.02) 


U 


0.70 








4780 


31.462 


0.351 


0.090(0.046) 


97.3 


3 


0.09(0.01) 


N 


0.82 


5.56t°i^ 






4781 


31.466 


0.185 


0.255(0.054) 


103.6 


3 


0.19(0.08) 


U 


0.70 








4794 


31.516 


0.449 


0.184(0.048) 


83.7 


3 


0.11(0.04) 


N 


0.91 


A Q4+0-42 






4811 


31.580 


0.227 


0.186(0.048) 


115.7 


1,3 


0.14(0.06) 


T 


0.57 


7 no+0''* 

'■'^^--0.62 


7.13(0.66) 


4814 


31.584 


0.205 


0.218(0.052) 


114.9 


3 


0.15(0.04) 


T 


0.61 


fi Sa+0.82 
O-88-0.56 


7.07(0.66) 


4826 


31.608 


0.171 


0.120(0.042) 


105.4 


3 


0.10(0.03) 


U 


0.60 







References. — 1: HCO + I IShirlev et al]|2013ft : 2: CS (Y. Shirley 2012, private communication); 3: NH3 I IDunliam et al.ll20lil 'l 

Note. — Table |3] is published in its entirety in a machine-readable format in the online journal. A portion is shown here for guidance 
regarding its form and content. 
'^ IRosolowskv et al.l 11201^') 

Flux density and uncertainty within a 40" aperture, corrected by the factor of 1..5 it 0.15 from lAguirre et al.1 ll201ir 
*^ N — near; F =^ far; T — tangent point; U — unconstrained distance 

Integrated posterior DPDF on the side of tZtan containing d^, . Larger values indicate higher certainty in the KDA resolution. 
^ Maximum-likelihood distance; the distance where the posterior DPDF is largest. Not listed for unconstrained sources. 

Weighted-average distance; the first moment of the posterior DPDF. Only listed for sources at the tangent point. 



pixels were removed from consideration (type 2). 
These additions to the automated selection criteria 
led to fewer sources (only 28%) requiring by-eye re- 
moval, primarily due to poor star-subtraction (type 
1) or complex emission structures that caused mor- 
phological mismatch (type 4). 

5.1.3. Source Properties 

The final source list contains 770 BGPS objects, 
and is presented in Table O EMAF-selected BGPS 
molecular cloud clumps are not drawn uniformly 
from the BGPS catalog. Comparisons of Galactic 
latitude and 40" flux density distributions between 
this sample and the full Bolocat (within the spatial 
limits defined above) are shown in FigurelH The lat- 
itude distribution of this sample follows that of the 
BGPS as a whole, including peaking below b = 0°. 
The offset is related to the Sun's ve rtical displace- 
ment abo ve the Galactic mi dplane (iSchuUer et al.l 
[2009: Ros olowskv et al.ll2010( ). The only significant 
deviation is near b = 0°, where locally bright 8- 
/xm emission along the midplane, excited by H II 
regions and OB stars, obscures more distant molec- 
ular cloud clumps. The BGPS 40" flux density his- 
tograms (Fig. |6^) show that this sample contains, 
on average, brighter sources (median = 0.252 Jy) 
than the full Bolocat (median = 0.135 Jy). There 
are two likely origins of this bias. First, sources 
must have a v^^ measurement from a dense-gas 
tracer; the HCO^ detection fraction, in particular, 
is a strong function of BGPS fiux density (< 20% for 



5"^, < 0.1 Jv: iShirlev et al.ll2013D . Second, the se- 
lection criteria excluded sources with very low con- 
trast or whose morphology does not correspond to 
dark regions in the GLIMPSE maps. Faint BGPS 
sources have low optical depth (5^ ^ = 0.1 Jy cor- 
responds to Tg w 0.07), and would be difficult to 
distinguish against the variable Galactic 8-/.tm back- 
ground. 

The distribution of measured mid-infrared source 
contrast is shown in Figure [Blc, and has a median of 
0.19. For images without the IRAC scattering cor- 
rection, this value corresponds to an uncorrected 
median contrast of 0.15 (see Appendix IBJ) . consid- 
erably l ower than the min i mum contrast (C w 0.20) 
used by iPeretto fc FuUeil ()2009[ ) in their catalog of 
SpitzerlKDCs (which did not correct for IRAC scat- 
tering in the same manner). The majority of our 
sample consists of "low-contrast" sources that are 
missing from published catalogs of Galactic IRDCs. 

5.2. Source KDA Resolutions 
5.2.1. Distance Estimates and Constraints 

We derived posterior DPDFs for the EMAF- 
selected BGPS sources by multiplying the kinematic 
distance DPDF by the two priors, and normaliz- 
ing to unit total probability. By design, the prior 
DPDFs are broader than the peaks in DPDFkin, so 
the resulting maximum-likelihood distances gener- 
ally do not differ from the simple kinematic dis- 
tances by more than ~ 0.1 kpc. To gauge the 
strength of the KDA resolution, two statistics were 
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Fig. 7. — Comparison of the fuU-wridth of the 68.3% error 
bar (FWes) against the integrated DPDF on the d-^^ side 
of dtan (^ml)- "^^b vertical dot-dashed line represents the 
empirical cutoff at Pj^^ = 0.78, and the horizontal dashed 
line marks FWes = 2.3 kpc. Objects shown in gray are below 
the PpjL cutoff and are more than 1 kpc from the tangent 
point. We defined FWes < 2.3 kpc as the criterion for a 
"well-constrained" distance estimate, which encompasses the 
objects in the bottom left corner, as well. 



defined: the maximum-likelihood probability (P^l) 
as the integrated posterior DPDF on the d^^^ side of 
the tangent point, and the full width of the 68.3% 
maximum-likelihood error bar (FWes). The ranges 
of these statistics are 0.5 < P^^^ ^1-0 and 0.2 kpc < 
FWgs ^ 15 kpc, and a comparison between them for 
each object is shown in Figure [T] The nature of a 
double-peaked DPDFkin leads to the sharp change 
in the distribution of FWgg near P^^^ =0.78 (verti- 
cal dashed line) . When the ratio of the peak prob- 
abilities of the kinematic distance peaks in the pos- 
terior DPDF becomes < 3, the error bars must in- 
clude both to enclose sufficient probability. For ob- 



jects with PjjL ^ 0.78, the maximum-likelihood er- 
ror bars enclose a single kinematic distance peak. 
We consider this set to have well-constrained dis- 
tance estimates, and note that FWes £ 2.3 kpc 
(horizontal dashed line). Objects with full- width 
error bars less than this value and are below the 
P„L cutoff (lower-left corner of Fig. [7]) are generally 
within ^ 1 kpc of the tangent distance. Because of 
the limited distance range available to these sources, 
their distance estimates should also be considered 
well-constrained. Combining these sets, we adopted 
FWes < 2.3 kpc as the criteria for well-constrained 
distance estimates. 

Objects within a kiloparsec of dtan have posterior 
DPDFs that are oftentimes asymmetric, and d^^ is 
not the best single-value representation of the dis- 
tance. For these objects, we assigned them to the 
tangent distance group, and used d (and associated 
uncertainty) in the analysis that follows. A total 
of 618 sources in this sample have well-constrained 
distance estimates (80%). The KDA resolutions for 
these objects are recorded as "N" (near), "F" (far), 
or "T" (tangent distance) in column 8 of Table [SJ 
d^L is listed in column 10. The remaining objects 
are recorded as "U" (unconstrained) and have no 
distance estimate listed. Foi^the tangent group, the 
weighted-average distance (d) is listed in column 11, 
and is the preferred distance representation for these 
objects (d„L is shown for comparison only). DPDFs 
for all sources are available in the BGPS archive. 

5.2.2. Heliocentric Distances 

For the remaining discussion, we consider only 
the 618 mid-infrared-dark BGPS sources whose dis- 
tances are well-constrained (as defined above). Of 
this set, 70 were placed beyond the tangent point, 
with another 25 near dtan, indicating the significant 
possibility of detecting molecular cloud clumps at 
the far kinematic distance using mid-infrared ab- 
sorption. The comparisons of latitude distribution, 
BGPS 40" flux densities, and mid-infrared contrast 
between the near and far subsets are shown in Fig- 
ure m For the purposes of this discussion, objects 



16 



Ellsworth-Bowers ct al. 



80 




Near — . 
Far + Tan — 


^60 




1 


1 


1-1 


■ 


T3 




: 






d 40 


■ 


! 


■ 


0) 
Q. 


p 


J 1 




^ 20 


■ 


JA ■ 







. . . .1. . . . . . . n-, . , 



100 



.E 80 

j3 



A 60 
en 

o 40 

Q. 



-1.0 -0.5 0.0 0.5 1.0 
Galactic Latitude [deg] 



20 




(b) 



m 



Near — 
Far + Tan — 



1 



^ 



0.1 1.0 

BGPS 40" Flux Density [Jy] 



60 

g 
i 50 

CO 

^ 40 

o 
o 

(M 30 

o 

Q. 

z 10 



(c) 



% 



it 



Near 
Far + Tan 



In 



. \n. 



0.0 0.2 0.4 0.6 

Mid-Infrared Contrast 



Fig. 8. — Comparison of source properties for objects with "near" vs. "far" KDA resolutions. Sources placed at the near 
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at the tangent distance are grouped with those at 
the far kinematic distance. The latitude distribu- 
tions (panel a) are very similar, with the near group 
being slightly wider, owing to the latitude-limiting 
effect of DPDFh2 • 

The histograms of BGPS 40" flux densities 
(Fig. [8j)) show that the far subset has a flatter dis- 
tribution with a higher median than the near set. 
Since the source list for this study is mid-infrared- 
contrast limited, we do not expect to see low flux- 
density BGPS sources at the far distance; the low 
column density would not produce enough attenu- 
ation to be seen behind the significant foreground 
emission. The expected contrast as a function of 
heliocentric distance (represented by /fore) may be 
computed by combining Equations ([7]), ©, and (|lip 
into 

C=(l-/forc)(l-e-^^iO , (12) 



2.0 



where /omaf and /„, 
lent to /niin and (L. 



from Equation ([7]) are equiva- 
J from Equation pTj) . respec- 



tively. A source with larger flux density may be at 
a farther dg and still meet the contrast selection 
criterion. 

Indeed, sources at the far kinematic distance have 
a lower median mid-infrared contrast (C = 0.11) 
than those at the near kinematic distance (C = 
0.22), as shown in Figure [He. Of the 313 objects 
with C > 0.2, only 12 (4%) are placed at the far 
kinematic distance, reinforcing the notion that dark 
mid-infrared absorption features must lie relatively 
nearby. For comparison, of the 63 sources with 
C < 0.1, 40 (63%) were placed at the far kinematic 
distance, indicating that the majority of EMAFs 
with very low contrast are at or beyond the tangent 
point. 

An interesting empirical predictor of KDA reso- 
lution is shown in Figure [9] The ratio of resolved 
heliocentric distance over dtan is plotted against the 
ratio of BGPS 40" flux density over the mid-infrared 
contrast. For BGPS data, there appears to be a 
boundary at S^ j/C ~ 2.5 that divides KDA reso- 
lutions. The exact value of this cutoff is dependent 
upon the (sub-)millimeter survey used, and there 
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Fig. 9.— KDA resolution versus the ratio of BGPS 40" 
flux density to mid-infrared contrast. The upper region rep- 
resents the far kinematic distance, and the lower the near; 
the gray shaded region illustrates the band around dtan- 
Black dots mark the sources with well-constrained distance 
estimates; the s ubset of cyan squares are sources within 
W43 (see TO. 1. 111 . The vertical dot-dashed line is drawn at 
S-^ j/C = 2.5. (A color version of this flgure is available in 
the online journal.) 



exists scatter across the boundary. It nevertheless 
suggests an additional means for KDA resolution 
when DPDFomaf fails to return a well-constrained 
estimate. 

5.2.3. Galactocentric Positions 

With well-constrained distance estimates, it is 
possible to construct a face-on view of the Milky 
Way. Sources with well-constrained KDA reso- 
lutions are plotted atop a reconstruction of the 
Milky Way from Spitzer data in Figure [TUlfR. Hurt: 
NASA/JPL- Caltcc h/SSC) using either d^^ or d as 
described in i i5.2.1l For clarity, the error bars, which 
account for small deviations from circular motion, 
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Fig. 10. — Face-on view of the Milky Way for sources 
with well-constrained KDA resolutions, plotted atop an 
artist's rendering of the Milky Way (R. Hurt: NASA/JPL- 
Caltech/SSC) viewed from the north Galactic pole. The im- 
age has been scaled to match the /?o used for calculating kine- 
matic distances. The outer dotted circle marks the Solar cir- 
cle, and the inner dotted circle the tangent point as a function 
of longitude. The dashed circle at Rgai = 4 kpc outlines the 
regio n influenced by the long Galactic bar l IBeniamin et al] 
I2005I V corresponding to the hashed regions in Fig. [5] The 
straight dashed gray line marks I = 30° as a guide. Various 
suggested Galactic features are labeled. For clarity, distance 
error bars are not shown. (A color version of this figure is 
available in the online journal.) 



are not shown. Some spiral structure is evident in 
the map of BGPS sources, notably portions of the 
Sagittarius Arm at £ > 35°, the Scutum-Centarus 
Arm / Molecular Ring at £ < 30°, and the local 
arm / Orion spur with i n abo ut a kiloparsec of the 
Sun (jChurchwell et al.l [2009f ). The kinematic re- 
strictions on our sample led to the absence of ob- 
jects within a « 3.5 kpc radius of the Galactic center 
(dashed circle in the figure). Face-on views of the 
Galaxy derived from kinematic distances will not 
show narrow spiral features (like those in the back- 
ground image) because of the local virial motions 
of individual molecular cloud clumps within larger 
complexes. Galactocentric positions are therefore 
"smeared-out" by approximately ±0.4 kpc about 
the true kinematic distance for the complex as a 
whole. Each dot in the figure, however, should be 
thought of in terms of its DPDF, where the kine- 
matic distance peaks have a FWHM of 1-2 kpc. 

KDA resolutions also allow the derivation of the 
vertical distribution of sources about the Galac- 
tic midplane. Vertical position is particularly af- 
fected by the KDA for higher- latitude sources {\h\> 
0?4). Calculation of vertical height (z) requires a 
proper accounting of the Sun 's ~ 25 pc vertical off- 
set above the Galactic plane (jHumphrevs fc LarsenI 
119951 iJuricet al.ll2008[) . The small scale height of 
Galactic molecular gas can lead to incorrect infer- 
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Fit centroid = +6.7 ± 0.9 pc 
FWHM = 49.8 + 1.6 pc 
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Fig. 11. — Vertical distribution of sources about the Galac- 
tic midplane. The filled gray histogram shows the distribu- 
tion, while the black line represents a Gaussian fit to the 
histogram. 
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Fig. 12. — Derived vertical position of sources versus he- 
liocentric distance. Diagonal cyan dashed lines represent 
the nominal |b| = 0?5 limit of the BGPS at f = 30° for a 
vertical Solar offset above the Galactic midplane of 25 pc. 
Horiz ontal dot-dashed lines mark the FWHM of the ^^CO 
layer IIBronfman et aLlll988l ). (A color version of this figure 
is available in the online journal.) 



ences about the vertical distribution of dense gas in 
the disk if z positions are calculated directly from 
Galactic coordinates without a correction for the So- 
lar offset. The matrix needed to transform {£,b,d) 
into (i?gai, </<, z) is derived in Appendix ICl 

The vertical distribution for the set of well- 
constrained BGPS sources is shown in Figure [TTJ 
A Gaussian fit to the distribution yields a half- 
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TABLE 4 
Maseh Sources for Distance Comparison 



Source 
Name 



(°) (°) (kms-i) 



Distance 
(kpc) 



N^ 



Ref. 



G23.0-0.4 


23.01 


-0.41 


81.5 


4 6+"'* 

^■°-0.3 


G23.4-0.2 


23.44 


-0.18 


97.6 


K q+1-4 


G23. 6-0.1 


23.66 


-0.13 


82.6 


•^■^-0.4 


W51 IRS2 


49.49 


-0.37 


56.4 


5.1+fl 


W51 Main'' 


49.49 


-0.39 


58.0 


c 4+0.3 



6 


1 


2 


1 


2 


2 


1 


3 


1 


4 



References. — (1) iBrunthaler et all (l2009fl: (2) 

iBartkiewicz et all l|2008f ): (3) IXu et al.l l|2009l ): (4) ISato et all 
||20TOD 

'^ Number of EMAF-selected BGPS sources within 15' and 
10 km s~^ of the maser location. See Fig. [T2]for the comparison. 
'^ H2O maser; all others are CH3OH masers 



width at half maximum of 25 pc, and a posi- 
tive centroid offset of 7 pc. Thi s scale height is 
appro ximately half that found by iBronfman et al.l 
(jl988[) for ^^CO. This narrow result may, however, 
be a result of the limited Galactic latitude cov- 
erage of the BGPS. Analysis of the recent com- 
pact source catal og from the A = 87 /xm AT- 
LAS GAL survey ([Gontreras et al.ll2013() . which ex- 
tends to |6| = 1°, shows that w 20% of their ob- 
jects lie outside the BGPS latitude limits. To gauge 
the effect of limited latitude coverage, the derived 
z are plotted against heliocentric distance in Fig- 
ure [12 with l&l = 0?5 shown for i = 30° (the 
limits rotate to slightly more positive z for larger 
£). For the region dg < 6 kpc (which contains 
more than 80% of this sample), the BGPS does 
not fully probe the FWHM of the ^^CO distribu- 
tion (dot-dashed lines). Other indicators of a larger 
scale height for star-for niation regions include ^^CO 
clouds from the GRS (jRoman-Duval et al.l [20091) . 
which have a FWHM w 8 pc, and Galactic H II 
regions (FWHM w 100 pc: IAnderson et al.|[20Tl) . 

5.3. Distance Comparisons with Other Studies 

The quality of KDA resolutions for EMAF- 
selected BGPS sources was characterized through a 
comparison of distance estimates with values from 
the literature. In particular, comparison sets were 
chosen that used mostly orthogonal methodologies 
so that distance comparisons are largely free of cor- 
related effects. The three sets described below are 
the use of maser parallax measurements, H I ab- 
sorption features associated with molecular clouds, 
and ncar-infrarcd extinction measurements. 

5.3.1. Maser Parallax Distances 

Maser parallax measurements towards regions 
of high-mass star formation provide absolute dis- 
tance validation comparisons. T he Bar and Spiral 
Struc ture Legacy Survey (BeSSeL: IBrunthaler et al.l 
120 111 ) is conducting ongoing VLBI parallax mea- 
surements of GH3OH and H2O maser emission in 
star-forming regions across the Galactic plane. Such 
measurements provide very accurate distances out 



to dg ^ 10 kpc, but the present overlap between 
published results and the BGPS is small (see Ta- 
ble [4| for the comparison set of maser sources used) . 

The comparison set was defined as objects from 
our sample whose angular separations and veloc- 
ity differences were < 15' and < 10 km s~^, re- 
spectively, from those of a published maser. These 
masers tend to be in regions of high-mass star for- 
mation, and such regions are on order 0?25 in size. 
The velocity limits are related to the spread of 
virial velocities within such regions. A collection 
of 12 BGPS objects were associated with one of five 
masers; the distribution is noted in Table [H To 
visualize the distance comparison, distances from 
Table [3| for each BGPS object are plotted against 
the measurements from the BeSSeL literature as 
magenta triangles in Figure 1131 The gray dashed 
lines in the left panel represent ±1 kpc error mar- 
gins, used for qualitative purposes. An object falling 
outside this region is said to have a "mismatching" 
distance estimate. For clarity in the figure, only 
mismatching objects have error bars shown; hori- 
zontal bars are from Table [31 and vertical bars are 
from the reference. The right panel is a zoom-in on 
the Ad = ±1 kpc region of the left panel {i.e. within 
the dashed lines). 

For the maser comparison set, only three sources 
fall outside the ±1 kpc region. Two BGPS ob- 
jects (overlapping triangles in Fig. [T3|) are associ- 
ated with the maser source G23.66-0.13, which has 
a parallax distance that d isagrees with the derived 
(near) kinematic distance. IBartkiewicz et ahl (l2008f ) 
find that this object has a proper motion consistent 
with the parallax distance and the assumption of 
a flat rotation curve, but has a « 35 km s~^ pe- 
culiar motion toward the Galactic center. This ra- 
dial streaming motion makes its kinematic distance 
appear larger, and provides a cautionary example 
of the effects of non-circular motion on kinematic 
distance methods. The other mismatching source 
has a DPDF distance estimate skewed away from 
the simple kinematic distance due to the sharply- 
peaked DPDFemaf causcd by its bright millimeter 
flux density (5^ ^ = 3.9 Jy). The W51 region lies 
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Fig. 13. — Left: Comparison of KDA resolutions derived from the DPDF witli publislied distance estimates. Gray dashed 
lines represent ±1 kpc away from equality. Maximum-likelihood horizontal error bars are shown for GRS- and maser-associated 
sources lying outside this region. Vertical error bars for maser-associated sources come from Tabled Sec text for discussion 
of systematic offset for NIREX sources. Right: A zoom-in on the region ±1 kpc from distance equality to better visualize the 
distance comparison for these objects. The vertical axis represents the comparison set distance minus the heliocentric distance 
from the DPDF. (A color version of this figure is available in the online journal.) 



near the tangent point, so correct DPDF distance 
placement for these objects merely implies that the 
region's circular velocity is consistent with the ro- 
tation curve. The remaining EMAF-selected BGPS 
objects in this set have KDA resolutions that agree 
with the trigonometric parallax distance. 

5.3.2. Galactic Ring Survey KDA Resolutions 

For a larger distance comparison set, we used 
the KDA resolutio ns from the BU-FCR AO Galactic 
Ring Survey (GRS: lJackson et al.l[2006l ). By match- 
ing ^'^CO(l-O) emission morp hology and spectra 
with H I absorption features, IRoman-Duval et al.l 
(J2009D estimated the distances to some 750 molec- 
ular clouds in the inner Galaxy. Those authors 
used a combination of H I self-absorption (HISA j3 
and 21-cm continuum absorption features to pos- 
itively resolve the KDA. These techniques exploit 
the spectroscopic dimension of H I surveys, where 
cold atomic hydrogen within dense molecular gas 
absorbs line emission from warm gas at the same 
Wlsr on the far side of the Galaxy or continuum ra- 
diation from H II regions. Distance resolutions from 
this method are subject to uncertainties from non- 
circular and radial streaming motions, but are di- 
rectly comparable with the KDA resolutions of the 
DPDF method. 

EMAF-selected BGPS objects were associated 
with cataloged ^^CO clouds based on spatial and 
kinematic proximity. The association volume was 

^'^ Absorption features caused by cold neutral hydro- 
gen within molecular clouds are also called "narrow" 
self-absorption (HINSA) to distinguish them from the 
broader self-absorptio n features of diffuse H I clouds (c./. 
ILi fc Goldsmithll2n03r) . 



defined as a circle of radius 10^ (approximately the 
median size of a GRS cloud; iRoman-Duval et al.l 
I2009D . and a velocity spread equal to the ^^CO 
velocity dispersion, centered on the (f, &,Ulsr) co- 
ordinates from the GRS catalog. A total of 213 
EMAFs lie within the association volume of one or 
more GRS clouds. To ensure the accuracy of the 
associated GRS KDA resolution, BGPS flux den- 
sity maps were compared to both ^''CO intensity 
maps integrated over the velocity of the appropri- 
ate dense-gas tracer from Table [1] and H I 21-cm 
"on"-"off" integrated intensity (HISA) maps. For 
a handful of sources (~ 7%), a strong HISA signa- 
ture was present within the BGPS source contour 
even though the associated GRS cloud was placed 
at the far kinematic distance. None of these objects 
was listed as having a 21-cm continuum source, so 
the absorption signature is the result of cold gas at 
the near distance. These discrepant objects may be 
the result of line-of-sight confusion, a slight veloc- 
ity ofl^set between the parent cloud and the BGPS 
object, or incorrect association with a ^^CO cloud. 
Whatever the cause, the KDA resolution for the as- 
sociated GRS cloud was amended to "near" to re- 
fl ect the HISA signatu re. 

iRoman-Duval et all used the iGlemensI (|1985l ) 
curve to derive heliocentric distances, and differ- 
ences in rotation curve definition can cause distance- 
comparison discrepancies unrelated to the KDA (see 
Fig. [S]). To eliminate potential systematic effects, 
the KDA resolution and v^sr of each associated 
GRS cloud wer e mapped to a ne w heliocentric dis- 
tance using the lReid et al.l (|2009l ) rotation curve. In 
the comparison between the GRS-derived distance 
and those from the DPDFs (black dots in Fig. [T^. 
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nearly 92% of our distance resolutions match those 
of the GRS. This success rate is robust for the entire 
EMAF set, as enforcing a minimum mid-infrared 
contrast of C > 0.15 only increases the matching 
rate to 94%. 

The 17 BGPS objects with mismatching distance 
resolutions are shown with horizontal error bars 
from the DPDFs. Those in the upper- left of Fig- 
ure [13] have a l arge apparent mid-infrared absorb- 
ing column, but iRoman-Duval et aD (|2009( ) did not 
find evidence of self-absorbing H I. Conversely, 
those in the bottom-right have HISA signatures 
but were placed beyond the tangent point by the 
posterior DPDF. Examination by eye of this latter 
group showed that the two sources farthest from 
the one-to-one line have slight underestimates of 
-'^MiR around the EMAF; the values in the postage- 
stamp image reflect dimmer nearby regions. The 
DPDFomaf in these cases selects the far kinematic 
distance peak despite the presence of HISA for these 
objects. 

There are four objects whose GRS distance esti- 
mate is 5.5 kpc > dg > 8 kpc and disagree with 
the DPDF-derived distance. These all lie within 
^1.5 kpc of the tangent point. Since the kinematic 
distance DPDFs do not have two fully distinct peaks 
in this region, the particular shape of DPDFgmaf 
can have a significant impact on the derived single- 
distance estimators. The mismatches are due to 
the source being near dtan, and not an incorrect 
KDA resolution. The remaining eleven mismatch- 
ing sources in the upper-left of Figure [T3] (GRS- 
far, DPDFcmaf-near) are moderately dark EMAFs 
(0.1 < C < 0.3) that show no signs of HISA at 
the velocity of the molecular cloud clump. About 
half of these lie at |6| > 0?4, and may not have 
enough H I backlighting at the far kinematic dis- 
tance for a H I SA si gnature to be visible; although 
iGibson et al.l (|2005D found self-absorption features 
out to more than |&| = 2° in the Canadian (H I) 
Galactic Plane Survey. For sources in this quadrant 
of the figure, it is unclear which kinematic distance 
is correct. The future application of additional prior 
DPDFs may solve the small number of conflicting 
KDA resolutions, but the present method achieves 
very good correspondence with other distance esti- 
mates for molecular cloud clumps. 

5.3.3. Near- Infrared Extinction Distances 

Using a technique for measuring three- 
dimension al near-infrare d Ga lac tic extinction 
NIREX; [Marshall et all [200l . [Marshall et aD 
2000) estimated the distances to over 1200 IRDCs 
identified by MSX in the inner Galactic plane 
(jSimon et al.ll200"6l hereafter S06). This approach 
compares the stellar colors of a section of sky with 
a Galactic stellar distribution model, and searches 
for sharp changes in color excess as a function of 
distance. Extinction measurements, like maser 
parallaxes, offer a kinematic-independent means of 
distance determination. 

MSX dark clouds have a typical size of about 



an arcminute, so BGPS sources lying within 60" 
of the ccntroid of a NIREX cloud were included in 
this comp arison set. While there are about 275 ob- 
jects from [Marshall et al.l (j2009( ) within the spatial 
bounds of this study, only 38 EMAF-sclectcd BGPS 
sources could be assoc iated with a NIREX cloud. 
iPeretto fc FJHIe^ (|2009l hereafter PF09) noted that 
only a quarter of MSX IRDCs appear in their cat- 
alog of Spitzer dark clouds for a variety of reasons. 
This selection effect, in combination with our re- 
quirement that an object be detected in one or 
more molecular line transitions, makes the number 
of matching clouds reasonable. 

Objects from the NIREX comparison set are 
shown as cyan diamonds in Figure |T3| Most of the 
points lie within 2 kpc of equality. Only one object 
has a wildly divergent distance estimate, G31.026- 



0.113 (BGPS #4653; < 



9.5 kpc), which 



has C = 0.4 and should have been detected with 
a strong near-infrared absorption signature at the 
near kinematic distance of d^ = 4.5 kpc. The col- 
lection of cyan diamonds with a systemic positive 
offset of 1.5 kpc warrants attention. There is a clus- 
ter of objects placed 2-4 kpc from the Sun. Most 
of these are at £ < 15°, and uncertainties in both 
the rotation curve and stellar model in that region 
may be contributing to the offset. Mismatching 
NIREX distances beyond rfg = 4 kpc have divergent 
dist ance estimates of or der the difference between 
the lClemeni (J1985I ) and lReid etai\ (|2009( ) rotation 
curves for objects at that velocity. 

6. DISCUSSION 
6.1. Kinematic Distance Discrimination 
6.1.1. EMAFs as Distance Discriminators 

The combination of millimeter- wave thermal dust 
emission observations with mid-infrared extinction 
is a powerful method for resolving the KDA for 
molecular cloud clumps. By starting from a catalog 
of (sub-)millimeter sources, this method is not lim- 
ited to mid-infrared-identified IRDCs (catalogs of 
which often have large minimum contrast). We are 
therefore able to include low-column nearby sources 
as well as more distant objects. The 93% success 
rate compared to distances resolutions by the GRS 
team indicates that EMAFs can provide a power- 
ful means for distance discrimination. Addition- 
ally, BGPS objects placed at the far kinematic dis- 
tance that agree with the GRS distance indicate 
that EMAFs are visible beyond the tangent point 
with sufficient backlighting. In comparison with 
the HISA KDA-resolution technique employed by 
IRoman-Duval et al.l (l2009l) . only 4% of BGPS ob- 
jects were placed at the near kinematic distance, 
yet had no evidence of a HISA signature. Applica- 
tion of additional prior DPDFs may help to resolve 
these discrepancies. 

The mid-infrared contrast distributions (Fig. \Ejp) 
of objects on either side of dtan clearly show that 
objects placed at or beyond the tangent point have 
lower collective contrast, and would likely not be 
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TABLE 5 
Effect of Dust Temperature on KDA Resolution 



Ti 


N ^ 


KDA Resolution*" 


GRS 


Comparison 


(K) 




N 


P 


T 


N 


Rate" 


15 


625 


416 


175 


34 


218 


77.5% 


20 


618 


523 


70 


25 


213 


91.9% 


25 


605 


547 


33 


25 


198 


91.9% 



" Number of well-constrained distance estimates 
'' N = near; F = far; T = tangent point 
'^ Distance matching success rate 



included in catalogs of IRDCs. These distributions 
are consistent with the notion that dark IRDCs 
(C > 0.2) are nearby. Since the matching rate 
between DPDF-derived distances and those of the 
GRS is nearly independent of mid-infrared contrast, 
the present method extends robust KDA resolution 
to EMAFs with lower contrast, roughly doubling 
the number of molecular cloud clumps for which 
well-constrained distances may be derived. 

In addition to improving upon the axiom "if 
IRDC then near" for KDA resolution, this method 
automatically accounts for the profile of the 8-/im 
intensity as a function of Galactic longitude (Fig.[2|). 
The morphological matching process does not con- 
sider dtan, and therefore offers a prior probability 
that is independent of the kinematic signature of 
a given object. The /fore liniit of visibility for a 
molecular cloud clump is simply a function of opti- 
cal depth (Equation fT^ : for instance, an object with 
S*! 1 = 0.3 Jy wiU have C > 0.05 for /fo^e < 0.76, but 
one with S'j^ ^ = 0.1 Jy will not meet this contrast 
threshold if /fore exceeds 0.33. 

Heightened star-formation activity, which pro- 
duces excess 8-^m emission in its immediate vicin- 
ity, does constitute a complicating factor in appli- 
cation of simple radiative transfer (Equation |B]). 
These regions strain the assumption of smooth, 
axisymmetric Galactic emission. As an example, 
we analyzed the distance resolutions of objects in 
the W43 region. W43 is defined here by 31?5 > 
e > 29?5, -0?5 > b > 0?3, and 80 km s'^ > 
f, a„ > 110 km s~^ (as in iNguven Luong et al.l 
120 lit) , and is marked by a white box in Figure [5) 
Of the 43 EMAF-selected BGPS sources with well- 
constrained distance estimates in this region, 9 are 
placed at or beyond the tangent distance by the 
DPDF method. This is a slightly higher rate than 
the general sample, but is not significant. The W43 
objects are plotted as cyan squares in Figure [9l 
and obey the empirical 5*^ ^/C = 2.5 limit for near 
versus far distance discrimination f §5.2.2|) . Only 8 
objects could be associated with a GRS-identified 
^■^CO cloud fi i5.3.2|) . and all but one have matching 
KDA resolutions; the outlier is one oHhe BGPS ob- 
jects near the tangent point, where d is used, caus- 
ing the > 1 kpc distance discrepancy. Although the 
observed Galactic plane consists of clumpy emission 
atop a more smooth Galactic emission pattern, the 
simple model used here still returns consistent KDA 
resolutions, even in more active regions. 



While the KDA resolutions for EMAF-selected 
BGPS objects compare favorably with previously- 
published distance estimates, they are still based 
upon the assumption of circular orbits about the 
Galactic center. As highlighted by the case of 
CII3OII maser G23.66-0.13, radial streaming mo- 
tions can have a significan t impact upon t he de- 
rived kinematic distances. lAnderson et al.l (|2012f ) 
presents a detailed analysis of uncertainties involved 
with the use of kinematic distances in the presence 
of non-circular motions. Future improvements in 
kinematic distance measurements will require a full 
three-dimensional vector model of Galactic motions. 

Throughout this analysis we used T^ = 20 K as 
the temperature for converting BGPS fiux densities 
into 8-fini optical depths. The effect of different 
assumed dust temperatures on KDA resolutions is 
not a priori predictable. Generating a new set of 
DPDFcmat based on different temperatures, how- 
ever, is straightforward. The results of KDA res- 
olutions from the posterior DPDFs and GRS dis- 
tance comparison statistics are shown in Table El fo r 
the range of Td found by iBattersbv et al.l (|2011| ) . 
A warmer dust temperature pushes more objects 
to the near kinematic distance to compensate for a 
smaller derived r^. Interestingly, the GRS success 
rate is unchanged for T^ = 25 K, but there are likely 
many unconsidered systematic effects at play. 

The use of EMAFs as prior DPDFs for kinematic 
distance discrimination is directly applicable to all 
current and future (sub-)millimeter surveys of the 
Galactic plane. The advantage of starting with 
a sample of continuum-identified molecular cloud 
clumps is that mid-infrared extinction may be used 
to resolve the KDA for many objects that may not 
be dark enough to be included in IRDC catalogs 
{e.g. S06, PF09). 

6.1.2. Use of Different Distance Estimates 

The distance probability density function 
(DPDF) formalism encodes all information about 
distance determinations for molecular cloud clumps, 
including likely distance and uncertainty. For some 
purposes, however, it is useful to have a single 
distance; §3.2! describes two possible options. In 
the use of the DPDFs produced here, it became 
apparent that different distance estimates were 
best applied to different situations. Examples of 
these situations are shown in Figure [14] to illustrate 
the difficulties encountered in extracting a single 
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Fig. 14. — Example DPDFs for tlirec possible cases. Shown 
in each panel are DPDFjjin (dashed gray), DPDFemaf (solid 
cyan), DPDFuj (dot-dashed gray), and posterior DPDF 
(solid black). The vertical dot-dashed line in each panel 
marks the tangent distance along that line of sight. The 
single-distance estimates are marked as black triangles (rfj^^ ) 
and magenta diamonds (d). Panel (a) represents a well- 
constrained KDA far from dtan; (b) shows a source near dtan; 
(c) illustrates a source with an unconstrained KDA resolu- 
tion. (A color version of this figure is available in the online 
journal.) 



distance from a DPDF. For each panel, the black 
triangle and magenta diamond mark d^^ and d, 
respectively. 

The most common situation has well-separated 
kinematic distance peaks {i.e. the molecular cloud 
clump is far from dtan) and the probability ratio of 
the two peaks is also large {i.e. P^^ — 0.78). The ex- 
ample BGPS #4484 (G030.629-00.029) is depicted 
in Fig. [Hb . These objects fall under the "well- 
constrained" condition described in ij5.2.11 whereby 
the maximum-likelihood error bars encompass only 
one kinematic distance peak. For this set of ob- 
jects, d„L is a reasonable collapse of the DPDF 
into a single value (with uncertainty). The second 
set of objects are those with kinematic distances 
within a kiloparsec of the tangent point. The kine- 
matic distance DPDF for these objects have a shal- 
low saddle feature at dtan, as seen in Fig. [Mb for 
BGPS #4357 (G030. 321-^00.292). Since the main 
peak of the posterior DPDF is not symmetric, d^^ 
is not a robust reflection of the distance estimate. 
Therefore, we recommend using d for these sources 
(listed in Table[3]for this class of object) to more ac- 
curately reflect the available distance information. 
So long as the full-width of the error bars is less 
than 2.2 kpc ( §5.2.ip . these sources are considered 
to have well-constrained distance estimates. The 
final set of sources are those not meeting either of 
the above criteria, such as BGPS #3352 (G024.533- 
00.182; Fig. [T4b). The kinematic distance options 
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Fig. 15. — Histograms of EMAF contrast as a function of 
IRDC identification. The black open histogram depicts ob- 
jects associated with the PF09 catalog; black dotted shows 
those associated with S06 (x5 for clarity); filled gray repre- 
sents objects associated with neither catalog. 

are well-separated, but the DPDFcmaf does not 
place a strong discriminatory constraint. If it is 
desirable to use the distances for these sources, we 
recommend Monte-Carlo sampling distances from 
the full DPDF in order to include all distance infor- 
mation about the source (see ^3.3^ . 

Regardless of the method used, however, care 
should be taken to properly propagate the uncer- 
tainty in the distance placement. If only sources in 
the first category are used, it is important to remem- 
ber that matching success rate with GRS-derived 
distances was sa 92%. This may be interpreted ei- 
ther as two sources in 25 were placed at the wrong 
distance or that there is a 92% confidence in each 
of the distance placements. 

6.2. EMAF vs. IRDC 

In this paper, we introduce the nomenclature 
"Eight-Micron Absorption Feature" (EMAF) for 
dust-continuum-identified molecular cloud clump 
whose emission morphology matches an absorption 
feature in mid-infrared maps of the Galactic plane. 
Many of these objects are quite dark (C > 0.2) and 
are identified in IRDC catalogs. To better under- 
stand the overlap between the EMAF and IRDC 
designations, we searched through both the S06 and 
PF09 catalogs to find the closest IRDC to the lo- 
cation of peak BGPS flux density. A total of 361 
(46%) EMAF-selected BGPS sources lay within the 
semi-major axis distance of the centroid of a cloud 
from one or both of the catalogs. 

By definition, IRDCs have a large mid-infrared 
contrast, but EMAFs are selected from thermal dust 
emission catalogs. As such, the two groups have 
different contrast distributions, as seen in the his- 
tograms of Figure [151 BGPS sources that are as- 
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sociated with an object in the PF09 or S06 cat- 
alogs are plotted as black solid and dotted lines, 
respectively. The filled gray histogram represents 
EMAFs not associated with any IRDC. The bulk 
of the non-IRDC objects have C < 0.2, once again 
confirming that the EMAF designation allows the 
use mid-infrared observations for KDA resolution 
for objects with low contrast. Of particular inter- 
est are the higher-contrast (C > 0.2) BGPS sources 
which do not appear in either IRDC catalog (filled 
gray). One striking example is the source G035.478- 
00.298 (BGPS #5631), shown in the lower-right cor- 
ner of Figure ^jp,. These objects suggest that it is 
easier to identify molecular cloud clumps from (sub- 
) millimeter data than to try to find intensity decre- 
ments in A = 8 fj,m images. 

The EMAF-derived KDA resolutions of the IRDC 
objects place only a small fraction at the far kine- 
matic distance (3% and 7% for S06 and PF09, re- 
spectively) . The measured EMAF contrast for such 
objects is generally C < 0.2, reinforcing the notion 
that dark IRDCs are nearby. Since the fractions of 
IRDCs placed beyond dtan are comparable to the 
GRS distance mismatch rate ( §5.3.2p . these subsets 
are not significant. 

6.3. Implications for Galactic Structure 
6.3.1. Galactic 8-jj,m Emission 

A quick glance at the GLIMPSE mosaics sug- 
gests that the Galactic distribution of 8-/im light 
cannot be described simply by a smooth, diffuse 
model. Distances derived using the assumption of 
smooth emission, however, compare favorably to 
those derived from 21-cm H I absorption (GRS; 
iRoman-Duval et al.ll20p9f) and near-infrared extinc- 
tion mapping (NIREX jMarshah et al.l2009( ). These 
results suggest that Galactic 8-/im emission may be 
primarily composed of a diffuse component punc- 
tuated by regions of active star formation. The 
R12 model of mid- and far- infrared emission is a 
greatly simplified reflection of the Galaxy, neglect- 
ing to account for individual small-scale features. 
Yet its match, spatially and spectrally, to existing 
Galactic plane observations indicates its power. The 
model, therefore, provides a solid basis for our dis- 
tance estimation technique based on using a broad 
distance prior to distinguish between kinematic dis- 
tance peaks. 

While the R12 model was constructed to broadly 
match the multi-band observations of the Milky 
Way, it is useful to also compare it with external 
galaxies. Spitzer obse rvations of SINGS galaxies 
(|Kennicutt et al.ir2003| ) yield a surface density of in- 
band emission. Figure [TBI depicts the IRAC Band 4 
surface brightness, normalized to Rgai = 8.5 kpc 
for a collection of 17 SINGS galaxies, with galaxies 
represented by their Hubble stage (TJ3- To enable 

^^ Hubble stage is a continuous numerical representa- 
tion of the Hubble type for a galaxy. T = corresponds 
to an SO galaxy, and the Milky Way (SBb-c) is T fa 4 
UBinnev fc Merriiield|[T998l . p.l55). 
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Fig. 16. — Azimuthally-averagcd surface brightness (A = 
8 ^in) as a function of galactocentric radius. A sample of 
17 SINGS galaxies is plotted in various line styles according 
to Hubble stage (T), and arc normalized to ijgai = 8.5 kpc. 
The thick black line depicts the surface brightness of the 
mid-infrared emission model (Fig. [l|. The depression in the 
model profile at low Rgai is due to the dust hole in the R12 
model ( ^4X3t . 

comparison, the R12 model was viewed externally 
(see Fig. [1]) and emission was annularly integrated, 
following the SINGS analysis; plotted as a thick 
black line. The slope of the Milky Way model at 
-Rgal ^ 4 kpc is comparable to the ensemble, but 
the drop in emission due to the dust hole carved 
out by R12 near the Galactic center docs not seem 
to match extragalactic observations. The model was 
designed to match observations from within the disk 
from the Sun's location, however, so discrepancies 
in the integrated 8-/im emission profile in the inner 
-Rgal < 3 kpc are likely not relevant. 

6.3.2. Spiral Structure 

Recent large (sub-)millimeter Galactic plane sur- 
veys are making it possible to trace the spiral struc- 
ture of the disk. Comparing the well-constrained 
KDA resolutions of BGPS sources with an artist's 
conception of the Galaxy based on Spitzer data 
(Fig, nop , glimmers of organization begin to appear. 
While the regions that trigger the collapse of molec- 
ular cloud clumps are likely very localized along spi- 
ral density waves, features in the face-on map of the 
Galaxy derived from kinematic distances are quite 
smeared out. Even small (« 5 km s~^) peculiar mo- 
tions can lead to « 0.4 kpc variations in heliocentric 
position, making it difficult to precisely trace the lo- 
cations of spiral features. 

Two major features suggest themselves from the 
data in Figure [TOl one near and one far. First, the 
nearby collection of sources at £ < 30° seem to form 
part of a round feature that extends into the fourth 
quadrant. This feature has been identified as both 
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the Molecular Ring (c./. iDame et afl 1200 ID and 
the S cutum-Ccntarus Arm (c.f. iDobbs fc BurkertI 
I2OI2D . It is not possible to distinguish between 
these postulates in the northern plane; careful dis- 
tance determinations for southern sources is re- 
quired. Mapping the exact location of where this 
feature meets the long Galacti c bar also depen ds on 
choice of rotation curve (the iClemensI 119851 curve 
places this collecti on of sources at the tangent dis- 
tance, whereas the lReid et al.l[2009l curve does not), 
and may also be influenced by non-circular motions. 
Parallax measurements to masers in this region will 
help establish a benchmark for the long Galactic 
bar, including position angle with respect to the Sun 
— Galactic Center line. 

The other feature to note in FigureHOlis the Sagit- 
tarius Arm beyond the tangent circle (smaller dot- 
ted circle) . These molecular cloud clumps are visible 
at the far kinematic distance because of backlighting 
provided by the Perseus Arm. By the same token, 
BGPS sources in the Perseus Arm are not visible 
as EMAFs due to both the large amount of 8-/.im 
light in the foreground and the lack of any signif- 
icant backlighting source. We note a collection of 
sources around £ ~ 30° which appear at the far 
kinematic distance. In the underlying image, there 
is a void between the Sagittarius Arm and the long 
Galactic bar. The two possibilities, therefore, are 
the existence of an arm structure at that location or 
that those sources were improperly assigned the far 
kinematic distance. The analysis of the W43 region, 
however, suggests that active regions do not signif- 
icantly misplace molecular cloud clumps at the far 
kinematic distance, so distinguishing between the 
possibilities is unclear. 

7. CONCLUSIONS 

We developed DPDFs as a new method for dis- 
tance determinations to molecular cloud clumps in 
the Galactic plane. Starting from a kinematic dis- 
tance derived from molecular line observations as 
the likelihood, prior DPDFs may be applied in a 
Bayesian manner to resolve the KDA. In this study, 
we used two external data sets as priors: mid- 
infrared absorption features, and the Galactic dis- 
tribution of molecular gas. 

The dust in molecular cloud clumps detected by 
(sub-) millimeter Galactic plane surveys should ab- 
sorb mid-infrared light and be visible against the 
broad diffuse PAH emission near A = 8 fim. Start- 
ing from the BGPS catalog of dust-continuum- 
identified molecular cloud clumps, we identified 
770 EMAFs in the Spitzer/ GLIMPSE mosaics. 
EMAFs may be thought of as generalized IRDCs, 
and are characterized by their selection from (sub- 
)millimeter data. With this collection of objects, 
simple radiative transfer arguments, and a model 
of Galactic mid-infrared stellar and dust emission, 
we developed a morphological matching scheme to 
compare dust emission and absorption. When us- 
ing the GLIMPSE mosaics to measure apparent ab- 
sorption features, it is imperative to account for 



scattering of light within the IRAC camera. This 
scattering, in concert with the instrumental calibra- 
tion method, means that diffuse emission will ap- 
pear brighter than it really is; bright emission will 
tend to fill in absorption features. The scattered 
light changes the apparent contrast of absorption 
features, in addition to any derived properties (such 
as optical depth, mass, etc.). 

Well-constrained KDA resolutions were obtained 
for 618 objects in this sample: 527 at the near kine- 
matic distance, 31 at the tangent point, and 70 at 
the far kinematic distance. To corroborate our dis- 
tance discriminations, we used VLBI niascr parallax 
measurements, KDA resolutions from the GRS, and 
near-infrared extinction distances as comparison 
sets. Of the 13 objects associated with maser paral- 
lax measurements, none had a discrepant KDA reso- 
lution. Distance comparisons with the GRS yielded 
a 93% success rate nearly independent of mid- 
infrared contrast. Comparison with the NIREX dis- 
tances showed only one discrepant KDA resolution, 
with the remainder being within identified system- 
atic effects. These comparisons illustrate the va- 
lidity of the present method, including the place- 
ment of some EMAFs at the far kinematic distance 
(approximately 12 distance-atched GRS sources are 
beyond dtan)- 

Approximately half of the set of EMAFs are as- 
sociated with an object from the IRDC catalogs of 
S06 and PF09. Objects associated with IRDCs are 
mostly relatively dark (C > 0.2); the remainder be- 
ing largely low-contrast (C < 0.2). Interestingly, 
there are a handful of moderately dark (C > 0.3) 
EMAFs that are absent from these IRDC cata- 
logs. This suggests that it is perhaps easier to iden- 
tify molecular cloud clumps first in (sub-)millimeter 
data, then investigate their mid-infrared properties. 

KDA resolutions for EMAFs from the BGPS cat- 
alog reveal hints of Galactic structure. Foremost, 
most detectable Galactic molecular cloud clumps 
are in the Molecular Ring / Scutum-Centarus Arm 
feature between the Sun and the Galactic center. 
The Sagittarius Arm outside £ = 30° is suggested 
by a collection of EMAFs beyond the tangent point, 
visible due to backlighting from the more-distant 
Perseus Arm. 

The derivation of DPDFs allows for probabilis- 
tic determination of distances to molecular cloud 
clumps across the Galactic plane. By introducing 
the concept of an EMAF, we were able to use the 
mid-infrared GLIMPSE data to resolve the KDA for 
many more sources than is possible with extant cat- 
alogs of IRDCs. Although this method applies only 
to - 10% of the BGPS catalog, the DPDF frame- 
work allows for the incorporation of additional prior 
DPDFs to expand the number of molecular cloud 
clumps with well-constrained distance estimates. 
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APPENDIX 

A. COMPUTING A BAND-AVERAGED DUST OPACITY: Spitzer IRAC BAND 4 

The apparent dust opacity of extinction features in broadband images is related to the dust opacity as a 
function of frequency and the spectrum of the light being absorbed. The IRAC Band 4 bandpass includes 
several distinct emission features from PAH molecules, as well as the complex behavior of the dust opacity 
near the 10-/im silicate feature. In this appendix, we derive a band-averaged dust opacity (k)^,^^^ for use 
with the GLIMPSE mosaics. 

Given the simple radiative transfer model of Equation ^, the intensity transmitted through a dust cloud 
is 



-'trans -'back ^ 



-i-h 



(Al) 



where /back is the background light (from the cloud to large heliocentric distance), and /trans is the transmitted 
light exiting the cloud on the near side {i.e. not including emission between the cloud and the observer). 
The band-averaged optical depth is related to the apparent dust opacity by (T)^^_^^^^ — S (K)ba„di where E is 
the mass surface density of dust. The band-averaged dust opacity is therefore 



= -T^ln 
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The ratio of intensities is computed from the band average over each quantity, 

-'trans \-'back,z.' ^ /band 
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where the v subscript denotes that quantity as a function of frequency, and /?ba„d (^) i^ ^^^ relative frequency 
response per unit power for the instrument. Since both averages are over the same response bandpass, the 

usual normalization terms cancel. 

The typical radius of an IRDC is small {^ 1 pc; iRathborne et al.l |2006[) compared to the accumulated 
path length (D) for the diffuse background (several kiloparsecs), so the intensity /back.i/ = J ji^ ds may be 
approximated by /back.i^ = jv D, where jj/ is the emission coefficient. Additionally, the relative response per 
unit power, /?band(^)j i^ proportional to {1/hv) S:^^^^(y), where Sy^_^_^^{v) is the relative response per photorP^. 
Canceling frequency-independent quantities and inserting the intensity ratio into Equation (jA2p yields the 
desired band-averaged dust opacity. 



{^).. 



In 



jv ' 5,„,,(i^).7, e-«- dv 



/-- 



•S'band [V) J^ dv 



(A4) 



For the emission spectra (j^), we used the dust emission models of iDraine fc Lil ()2007D . which contain a 
mixture of grain sizes in addition to a variable PAH mass fraction ((7pah)- The various emission spectra were 
derived by irradiating the dust with starlight intensity fields having a tunable minimum value t^min relative 
to the local interstellar radiation field {U = 1). 

The choice of dust opacity model (k^) has a nontrivial effect on the derived band-averaged opacity. Three 
different models were analyzed, and are shown in Table |6l First is the 0H5 model (dust grains with thin ice 
mantles, coagulating at 10^ cm~'^ for 10^ years; lOssenkopf fc Henninelll99^ used for determining the dust 
opacity at A = 1.1 mm for the BGPS. The remaining models were presented in lWeingartner fc Draind ()2001l 



' Obtained from http : //Irsa . ipac . caltech . edu/data/SPITZER/docs/lrac/callbratlonf iles/spectralresponse. 
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TABLE 6 
Computed {K)^^^^^ for IRAC Band 4 



Dust Emission^ (k), , 
9PAH ^min OH5'= WDOl-3.1'^ WDOl-S.S" 
(%) (cm2 g-l) 

0.47 0.1 1175 858 909 

1.0 1177 861 911 

10.0 1178 864 915 

2.50 0.1 1157 796 843 

1.0 1158 800 847 

10.0 1160 804 852 

4.58 0.1 1150 777 823 

1.0 1151 779 825 

10.0 1152 782 828 

^ Dust emission model from lDraine fc Lil l|2007f ) 
^ Extinction from lOssenkopf fc Hennind (|1994 Table 
1, Column 5) 

'^ Rv = 3.1 extin ction from WD Ol, with updated nor- 
malizations from IDraind (|2003f ) 

'^ Rv = 5.5, Case A, extinc tion from WD Ol, with 
updated normalizations from IDraind (|2003f ) 

hereafter WDOl), with updated normalizations given by iDraini ()2003rF1 . The second of the three models is 
the Rv = 3.1 Milky Way model, tried even though this value of the color excess per magnitude extinction is 
consistent with the diffuse ISM and does not hold for regions of dense gas. The final model utilizes case A 
for Ry = 5.5 (consistent with observations of molecular clouds), which sought to minimize the extinction 
differences between observation and model, while also including a penalty term to keep dust grain volume 
from exceeding abundance/depletion limits (WDOl). For each dust model, (k)^^^^^ was computed for three 
values each of Qp^g and f/min (Table [6]). 

The minimum value of the starlight intensity field has very little effect on the band-averaged dust opacity, 
meaning that the derived (K)-^^^^^ are valid for a wide range of environments. An order of magnitude change 
in the assumed PAH mass fraction cau ses only a 2% diff erence in the derived opacity for the 0H5 model, 
but the spread is 10% for the others. iDraine fc Lil (j2007f ) cite (/p^jj = 4.58% as best matching observations 
of the Milky Way. 

The 0H5 model is the preferred description of dust at (sub-)millimeter wavelengths (c./. iRathborne et al.l 
l2006t [SciiuUer et al.ll2009HAguirre et al.ll2011[). and has been used by p revious studies for estimating {K)^^^^^ 
for IRAC Band 4 images (jButler fc TanI 120091 : iBattersbv et al.ll2010f ). That model was computed from 
theory for coagulated grains (aggregates of smaller particles, with some voids) surrounded by an ice mantle. 
In contrast, the WDOl dust models utilized simple geometry (PAH molecules for very small grains, and 
graphite and olivine spheres for larger grains) and sought to fit a dust size distribution to parameterized 
observed extinction. 

The connection to obser ved extinctio n in th e infrared led us to choose the WDOl Rv = 5.5 model 
for this work. Following IDraine fc Lil ()2007l ). we used q^j^^ = 4.58% and Umin = 1.0 to compute 
Kg = (K)ba„d ~ ^25 cm^ g~^ of dust. For comparison, the corresponding value from the 0H5 model yields 
Kg = 1167 cm'^ g"-'^, a w 40% difference. We note that the preferred WDOl model predicts a BGPS dust 
opacity k-^ ^ — 0.272 cm^ g~^ of dust, approximately one quarter the value from 0H5. 

B. Spitzer IRAC SCATTERING CORRECTION FACTORS 

The IRAC camera on the Spitzer Space Telescope suffers from internal scattering within the detector 
arrays, particularly Bands 3 and 4. The s cattering is such t hat a fraction of the incident light on a pixel 
is distributed throughout the entire array (jReach et al.ll2005t IRAC Instrument Handbook). Image frames 
are converted into physical units (MJy sr~^) using point-source calibration data; point-source aperture 
photometry is therefore accurate because the calibration takes into account the light scattered out of the 
aperture and into blank sky pixels. Observed extended emission, however, has light from other areas of the 
array scattered into each pixel as well, and so will appear brighter than it really is given the point-source 
calibration. For the broad, diffuse emission of the Galactic plane in IRAC Band 4, there is a roughly constant 
positive offset of the measured intensity in each frame. 

For absorption features in the Band 4 images (EMAFs), however, the scattering cannot be corrected for 
by simple multiplicative aperture corrections. Because bright emission from surrounding regions is scattered 
into an EMAF, it will have a lower apparent contrast. To correct for this effect, an estimate of the scattered 
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light in a frame must be subtracted from each frame (S. Carey 2010, private communication), as was done 
in this study for a pixel-by-pixel comparison between GLIMPSE and synthetic S-fim images. 

Q uantities such as con trast and optical depth for EMAFs may be derived from the GLIMPSE mosaic s 
(e.g. iButler fc Tanll2009t PF09), but a scattering correction must be applied (c./. iBattersbv et al.l[201Cl( ). 
Because careful subtraction of scattered light is not always necessary for a given application, correction 
factors may be derived for quantities measured directly from the IRAC Band 4 data. In this appendix, we 
derive correction factors for mid-infrared contrast (C), 8-/im optical depth (r^), and foreground fraction of 
8-/im emission (/fore)- 

For regions of broad diffuse emission punctuated by dark clouds, the observed intensities /q and /i of the 
background and EMAF, respectively, are related to the actual intensities Sq and ^i by 

Io = So+X , and (Bl) 

h = S,+X, (B2) 

where X = ^ 5'o is the amount of scattered light, approximated by the fraction ^ = (1 — 0.737) = 0.263 of 
incid ent diffuse light sca ttered throughout the array, and 0.737 is the infinite-aperture correction for Band 4 
from lReach et alT(j2005D . Rearranging to compute the true intensities from observed quantities yields 

^o^^YTe) ' ^"""^ (B3) 

^1=^1-^0^^^. (B4) 

While the subtractive correction (Equation [BT| for the diffuse background is equivalent to a multiplicative 
correction (Equation IB3|) . correcting for the intensity within the EMAF is more complicated. To compute 
the true contrast of an EMAF, we begin with Equation ([TT|) : 

C.,„. = 1 - 1^ • (B5) 

Inserting Equations (|B3P and (jB4p . the true contrast becomes 

c.,.^(i-^),i + f) 

= C_(l + 0, (B6) 

where C^„„ is the quantity measured directly from the GLIMPSE images. The measured contrast will be 
smaller than reality, leading to an underestimation of optical depth and other quantities. 

With a measured EMAF contrast, it is possible to estimate the optical depth of a cloud or the foreground 
fraction of diffuse emission, given an assumption about the other. Equation (|12p may be rearranged to solve 
for either quantity in terms of the other. The optical depth (from which follows surface mass density) is 
given by 
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assuming a model that yields /fore (as in IButler fc Tan|[2009() . The true value of r^ will be larger by up to a 
factor of two for C^^^^ < 0.5. If, instead, the foreground fraction is desired given an external estimate of r^ 
(such as from (sub-)millimeter thermal dust continuum data; as in PF09), the true value is given by 

, 1 g.„o..(l+0 

/fore, true — J- -, _^ 

1 — e ^8 

i--^^^^^)(i + o-e 

1 — e a J 
= /toro,mcas (1 + C) ^ C • (B8) 

The actual foreground fraction will be smaller than that measured directly from GLIMPSE images, with the 
difference becoming less at large /fore- Any /fore.mcas ^ 0.2 maps to zero true foreground fraction, as negative 
values are not physical; such values arise from uncertainty in C and the derivation of r^ from (sub-) millimeter 
data. 
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C. THE VERTICAL SOLAR OFFSET AND CONVERTING (i^b^d^) INTO (i?gai,</>,z) 

Deriving Galactocentric positions of objects in the Milky Way requires a coordinate transformation of 
the triad {£,b,d^), where d^ is the heliocentric distance along the line of sight toward {£,b). The Galactic 
coordinate system was defined assuming the Sun is at the midplane of the disk (jBlaauw et al.lll960D. but more 
recen t studies have me asured a vertical solar offset of « 25 pc above the midplane (jHumphrevs fc LarsenI 
119951 iJuric et ani2008f). Since the ve rtical scale height of the molecular gas layer in the disk is small 
(HWHM Rs 60 pc: lBronfman et al]|1988[ ). neglecting to account for the solar offset may introduce a systematic 
bias in the derived vertical distributions of components of the Galactic disk. 

The coordinate transformation is done in cartesian coordinates. First the triad (£, ^jfig) is converted into 
local cartesian coordinates, where the x-axis is directed along the Sun — Galactic Center line, and z points 
north out of the plane, 

(xi \ / dg cos I cos b \ 
2/1 = dg siiil cos 6 . (CI) 

zi J \ d^., sin b J 

The local coordinates are transformed to the Galactocentric frame by (1) rotation by 180° in the x — y plane 
to place the -l-x-axis pointing away from the GC, (2) translation of the coordinate axes to place the origin 
at the GC, and finally (3) rotation by the angle in the x — z plane to place the -l-x-axis along the Galactic 
midplane (rather than along the Sun — Galactic Center line). The transformation may be written as 




cos 6* — sin 9 

10 

sine cose* 

1 



1 i?o 

10 

10 

1 



-10 
0-100 
10 
1 




(C2) 



where the lateral translation requires an augmented (affine translation) matrix. The resulting Galactocentric 
cartesian coordinates are 




i?o cos^- dg (cos/ cos 6 cosff + sinfo smO) 

— dg sin/ cos 6 
Rq sin 6' — dg (cos/ cos 6 sin 6' — sin 6 cos 61) 



(C3) 



The rotation angle 9 = sin~ {zq/Rq), where zq = 25 pc, corrects for the Sun's vertical displacement above 
the midpl ane. Galactocentric positions in the cylindrical coordinates (i?gai,0, z) may be extracted from 
Eqn. (jCSp in the usual manner. The rotation by 9 is most important for the derived Zgai, and has negligible 
effect on Rgai and (f>. 
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